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Abstract 

In this paper, we prove the existence of smooth initial data for the 2D free boundary 
incompressible Euler equations (also known for some particular scenarios as the water 
r"| wave problem), for which the smoothness of the interface breaks down in finite time into 

a splash singularity or a splat singularity. 
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I Introduction 

I.A Statement of the Problem 

In this paper, we prove that water waves in two space dimensions can form a singularity in 
finite time by either of two simple, natural scenarios, which we call a "splash" and a "splat" . 

The water wave equations (or 2D incompressible free boundary Euler equations) describe 
a system consisting of a water region il(t) C and a vacuum region \ Q{t), evolving as 
a function of time t, and separated by a smooth interface 

dVl{t) = {z{a,t) : a G M}. 
We write n^{t) = \ Qj^^j^ f]2(^^ ^ rpj^g g^i^ velocity v{x,y,t) G and the 

pressure p{x, y,t) eM. are defined for (x, y) G i}.{t). The fluid is assumed to be incompressible 
and irrotational 

V-v = 0, curl?; = mn{t), (LI) 
and to satisfy the 2D Euler equation 

[dt + {v-V^)]v{x,y,t) = -Vp{x,y,t)-{0,g) mn{t), (1.2) 

where > is a constant, and the term (0, g) takes gravity into account. 
Neglecting surface tension, we assume that the pressure satisfies 

P = P*{'t) at dO,{t), where p*{t) is a function of t alone. (1.3) 

Finally, we assume that the interface moves with the fluid, i.e., 

dtz{a, t) = v{z{a, t),t) + c*(a, t)daz{a, t), (1.4) 

where c^{a.^t) is an arbitrary smooth function of a,t (the choice of affects only the 
parametrization of d^{t)) and z{a,t) = (2;i(a, t), ^2(0;, t)). 

At an initial time t^, we specify the fluid region r2(to) and the velocity v{x, y, to) {{x, y) G 
f2(to)), subject to the constraint (I.l). We then solve equations (1. 1-1.4) with the given initial 
conditions, and we ask whether a singularity can form in finite time from an initially smooth 
velocity i'(-,io) and fluid interface dVt{to). 

The water wave problem comes in three flavors: 

• Asymptotically Flat: We may demand that z{a, t) — (a, 0) — )■ as a — )• ±00. 
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• Periodic: We may instead demand that z{a,t) — (a, 0) is a 27r-periodic function of a. 

• Compact: Finally, we may demand that z{a,t) is a 27r-periodic function of a. 

To obtain physically meaningful solutions in the Asymptotically Flat and Periodic flavors, 
we demand that 



where we regard 0,{t) as a subset of T x M, T = M/27rZ, in the Periodic case. 

In this paper, we restrict attention to periodic water waves, although our arguments can 
be easily modified to apply to the other flavors. (See Remark 1.5 below). 

Let us summarize some of the previous work on water waves. We discuss the real-analytic 
case later in this introduction. The existence and Sobolev regularity of water waves for short 
time is due to S. Wu [30]. Her proof applies to smooth interfaces that need not be graphs of 
functions, but [30] assumes the arc-chord condition 



The constant cac > is called the arc-chord constant, which may vary with time. 

The issue of long-time existence has been treated in Alvarez-Lannes [3], where well- 
posedness over large time scales is shown, and several asymptotic regimes are justified. By 
taking advantage of the dispersive properties of the water-wave system, Wu [32] proved ex- 
ponentially large time of existence for small initial data. 

In three space dimensions, Wu [31] proved short-time existence; and Germain et al [19], 
[20] and Wu [33] proved existence for all time in the case of small initial data. 

There arc several important variants of the water wave problem. One can drop the 
assumption that the fluid is irrotational. See Christodoulou-Lindblad [12], Lindblad [23], 
Coutand-Shkoller [16], Shatah-Zeng [28], Zhang-Zhang [35]. Lannes [21] considered the case 
in which water is moving over a fixed bottom. Ambrose-Masmoudi [4] considered the case 
where the equations include surface tension, and the limit where the coefficient of surface 
tension tends to zero. Lannes [22] discussed the problem of two fluids separated by an 
interface with small non-zero surface tension. Alazard et al. [1] took advantage of the 
dispersive properties of the equations to lower the regularity of the initial data. 

See also the papers of Cordoba et al. [13] and Alazard- Metivier [2]. 

In the case of large data for the two-dimensional problem (1. 1-1.4), Castro et al. in [10], 
[9] showed that there exist initial data for which the interface is the graph of a function, but 
after a finite time the water wave "turns over" and the interface is no longer a graph. For 
previous numerical simulations showing this turning phenomenon, see Baker et al. [5] and 
Beale et al. [6]. 

Next, we describe a singularity that can form in water waves. We start by presenting 
what we believe based on numerical simulations; then, we explain what we can prove. 



p{x,y,t) + gy = 0{l) in ^(t) 



and that 




z{a,t) - z{l3,t)\ > CAcW - ah a, /3 G R. 
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t=t, 



(a) The initial water region f2(to). 



(b) The water region ^l{ti) at a later time ti. 




(c) A "splash" forms at time t2 > ti. 



Figure 1: Evolution of a "splash" singularity. 



Our simulations show an initially smooth water wave, for which the fluid interface is a 
graph as in Figure 1(a). At a later time ti, the water wave has "turned over" as described 
in [10], [9], i.e., the interface is no longer a graph. Finally, in Figure 1(c), the fluid interface 
self- intersects at a single point ^, but is otherwise smooth. We call this scenario a "splash", 
and we call the single point at which the interface self-intersects, the "splash point". Beyond 
the time t2 pictured in Figure 1(c), there is no physically meaningful solution of (1. 1-1.4). 

Note that the arc-chord condition holds for times t < t2, but the arc-chord constant tends 
to zero as t tends to t2- 

The numerics that led us to Figures 1(a), 1(b) and 1(c) were performed using the method 
of Beale-Hou-Lowengrub [7] , with special modifications to maintain accuracy up to the splash. 
In this paper, we use the numerics only as motivation for conjectures, so we omit a detailed 
discussion of the algorithms used. Actual results from our simulations are shown in Figures 
3, 4 and 5. Figures 1 and 2 are cartoons. 

Now let us explain what we can prove regarding the splash scenario. Recall that [10], [9] 
already proved that a water wave may start as in Figure 1(a) and later evolve to look like 
Figure 1(b). In this paper, we prove that a water wave may start as in Figure 1(b), and later 
form a splash, as in Figure 1(c). 

We would like to prove that an initially smooth water wave may start as in Figure 1(a), 
then turn over as in Figure 1(b), and finally produce a splash as in Figure 1(c). To do so, 
our plan is to use interval arithmetic [24] to produce a rigorous computer- assisted proof that, 
close to the approximate solution arising from our numerics, there exists an exact solution 
of (1. 1-1.4) that ends in a splash. The stability result announced in [8, Theorem 4.1] is a 

^Here, we regard the fluid interface as sitting inside T x R; recall that our water waves are 27r-periodic 
under horizontal translation. 
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first step in this direction. We are grateful to R. de la Llave for introducing us to interval 
arithmetic and demonstrating its power. 

A variant of the splash singularity is shown in Figures 2(a) and 2(b). 




p = 




t = t, 



t = t, > t] 



(a) The initial water region 



(b) The "splat". 



Figure 2: Evolution of a "splat" singularity. 



The water wave starts out smooth, as in Figure 2(a), although the interface is not a graph. 
At a later time, the interface self-intersects along an arc, but is otherwise smooth. Again, no 
physically meaningful solution of (1. 1-1.4) exists after the time depicted in Figure 2(b). We 
call this scenario a "splat". In this paper, we prove that water waves can form a splat. 

The stability theorem announced in [8] shows that a sufficiently small perturbation of 
initial conditions that lead to the splash will again lead to a splash. We expect that the 
analogous statement for a splat is not true. 

We make no claim that the splash and the splat are the only singularities that can arise 
in solutions of the water wave equation. 

I.B Elementary Potential Theory 

To formulate precisely our main results, and to explain some ideas from their proofs, we 
recall some elementary potential theory for irrotational divergence-free vector fields v{x,y,t) 
defined on a region Q(i) C M? with a smooth periodic boundary {z{a, t) : a € M} for fixed t. 
We assume that v is smooth up to the boundary and 27r-periodic with respect to horizontal 
translations. We suppose that v has finite energy. 

Such a velocity field v may be represented in several ways: 

• We may write v = ^(p for a velocity potential 0(x, y, t) defined on il(t) and smooth up 
to the boundary. 

• We may also write v = V'^ip = {—dyip, dxip) for a stream function -0, defined on r2(t) 
and smooth up to the boundary. 

• The normal component of v at the boundary, given by 



Unormal{a,t) = v{z{a,t),t) ■ 




daz{a,t)\ 
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uniquely specifies v on Cl{t). Here, u"*" = {—U2,ui) for u = {ui,U2) G M^, and we 
always orient dQ{t) so that the normal vector {daz{a,t))^ points into the vacuum 
region \Q(t), 

The function Unormal{'^,t) satisfies 

/ Unormalioi,t)\daZ{a,t)\da = 0, 

but is otherwise arbitrary. 

Note that, because v has finite energy, (/) and tp are 27r-periodic with respect to horizontal 
translations. (Without the assumption of finite energy, (p and ip could be "periodic plus 
linear"). The functions (p and ^p are conjugate harmonic functions. 

• There is another way to specify v, namely 

,{x, y, t) = ^PV [ (^-^i(^.0.Z/-^2(^,t))| ^^^^ ^^^^^ ^^^^ ^ ^^^y 
27r Jt^\{x- zi{/3,t),y- Z2{l3,t))\^ 

for a 27r-periodic function co{P,t) called the "vorticity amplitude". See [5]. 

Formula (1.5) holds only in the interior of Q{t). Taking the limit as {x, y) — )• (zi{a, t), Z2{a, t)) G 
dO,{t) from the interior, we find that 

v{z{a,t),t) = BR{z,u;){a,t) + (1-6) 
where BR denotes the Birkhoff-Rott integral 

^ ' ' ^ 27r J^\{zi{a,t)-zi{P,t),Z2{a,t)-Z2{P,tW 

To see that v may be represented as in (1.5), (1.6), one applies the Biot-Savart law to a 
discontinuous extension of v from its initial domain Q,[t) to all of M^; to make the extension, 
one solves a Neumann problem in \ 0(t). 

Thus, our velocity field v admits multiple descriptions. Note that the description in terms 
of u is significantly different from the descriptions in terms of cp, ip and UnormaU because we 
bring in the Neumann problem on R^\i7(t) to justify (1.5) and (1.6). When dQ,(t) is a "splash 
curve" as in Figure 1(c), there is no problem defining cp and it is smooth up to the boundary, 
except that it can take two different values at the splash point, for obvious reasons. The 
same is true of ^p. Similarly, Unormaii^^^t) continues to behave well. 

However, there is no reason to believe that (jj{a^t) will be well-defined and smooth for a 
splash curve, since \ 0(t) is a somewhat pathological domain. Our numerics suggest that 
maxa |a;(a,t)| ~ where ts is the time of the splash. 

Let us apply the above potential theory to the water wave problem. A standard formula- 
tion of the problem [5] takes z[a, t) and tj(a, t) as unknowns. This has the advantage that at 
least we know where our unknown functions are supposed to be defined, which is more than 
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we can say for (f), ip and u. Standard computations (see e.g. [13, Section 2]) show that the 
water wave problem is equivalent to the following equations 



dtz{a, t) = BR{z, uj){a, t) + c{a, t)daz(a, t) 



(1.8) 



and 



dtio{a,t) = — 2daz{a,t) ■ dtBR{z,oj){a,t) 



- da ( TT^ — 12 ) (a> t) + da (c(a, t)uj{a, t)) 




+ 2c{a,t)daz{a,t) ■ daBR{z,u){a,t) — 2gdaZ2{oi,t). 



(1.9) 



Here, c(a, i) is a function that we may pick arbitrarily, since it influences only the 
parametrization of dfl{t). For future reference, we write down several standard equations 
that follow from (1. 1-1.4) by routine computation and elementary potential theory. 

Ax(l>{x,y,t) = AxtpiXjUji) = in Q{t); (f) and ip are harmonic conjugates. 



and n is the outward-pointing unit normal to dCl(t). 
il){x + 27r, y, t) = tp(x, y, t) and ^{x + 27r, y, t) = 4>{x, y, t) 
ip{x,y,t) = 0(1) as y ^ — oo 

V = V-^V in ^{t) 
dtz{a, t) = v{z{a, t),t) + c{a, t)daz{a, t) 



We may write u(a,t) to denote v{z(a,t),t). 

I. C Main Results 

Our main result is the following theorem. For the definition of a splash curve see Definition 

II. 1 in Section II. The interface shown in Figure 1(c) is an example of a splash curve. 

Theorem I.l Let z^{a) be a splash curve, where the splash point is given by z^{ai) = 2;° (02), 
ai ^ a2- Let Unormali^) ^ scalar function in H^{T), satisfying 




dMa,t) = ^ 



viz{a, i), t)p + c(a, t)v{z{a, t), t) ■ daz{a, t) - gy{a, t) + p*{t). (I.IO) 




(1.11) 



and 



'^lormal{(^-i.)^'^normal{(^2) < 0. 



(1.12) 
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Then there exist a time T > 0; a time-varying domain Q{t) defined for t G [0,T] and a 
velocity field v{x,y,t) defined for {x,y) G ^{t), t G [0,T] such that the following hold: 



il{t) and v{x,y,t) solve the water wave equations (I.I-I.4) for all t G [0, T]. 



(1.13) 



dft{t) is given as a parametrized curve {z{a,t) : a G M} 
with z{a, t) — {a, 0) 2'iT-periodic in a for fixed t. 



(1.14) 



z{a,t) - (a,0) G C([0, T], iJ^(T)) and v{z{a,t),t) G C([0, T], i7^(T)) 



(1.15) 



z{a, 0) 



Z (a) and Unormalia, 0) = Unormal{(^) M all a G M. 



(1.16) 



For each t G [0,r], the curve d^},{t) satisfies the arc-chord condition, 
but the arc-chord constant tends to zero as t 0. 



(1.17) 



This result was announced in [8]. 

To prove that "splash singularities" can form, we note that the water wave equations are 
invariant under time reversal. Therefore, it is enough to exhibit a solution of the water wave 
equations that starts as a splash at time zero, but satisfies the arc-chord condition for each 
small positive time. Theorem I.l provides such solutions. 

Since the curve touches itself it is not clear if the vorticity amplitude is well defined, 
although the velocity potential remains nonsingular. In order to get around this issue we 
will apply a transformation from the original coordinates to new ones which we will denote 
with a tilde. The purpose of this transformation is to be able to deal with the failure of 
the arc-chord condition. Let us consider the scenario in the periodic setting and then the 
transformation defined by z{a, t) = P{z{a, t)) where P is a conformal map that will be given 
as: 



and the branch of the root will be taken in such a way that it separates the self-intersecting 
points of the interface. We will also need that the interface passes below the points (ibTr, 0) 
(or, equivalently, that those points belong to the vacuum region) in order for the tilde region 
to lie inside a closed curve and the vacuum region to lie on the outer part. See Figures 3 
and 4. Here P{z) will refer to a 2 dimensional vector whose components are the real and 
imaginary parts of P{zi -\- iz2)- Its inverse is given by 



In this setting, P ^(z) will be well defined modulo multiples of 27r. 

Remark 1.2 Note that P{z) is periodic such that P{z + 2A;7r) = P{z). Moreover, P{z) is 
one-to-one in the water region and single-valued except at the splash point. 




P ^{w) = Hog 



1 -f- iw^ 



) 



= 2arctan('«;^) for w eC 
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0.4 
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Figure 3: Splash singularity at times t = (Red - splash), t = 4 ■ W ^ (Blue - turning) and 
t = 7- 10^3 (Black - graph). 

Remark 1.3 Although the transformation to the tilde domain is convenient, the real reason 
for Theorem I.l is that the potential theory inside the water region does not go bad as we 
approach the splash even though it goes bad in the vacuum region. 

We define the following quantities: 

ip{x, y, t) = il}{P~^{x, y),t), cl>{x, y, t) = (p{P~'^{x, y), t), v{x, y, t) = V(j){x, y, t), 

$(a, t) = (t>{z{a, t),t), §(a, t) = ip{z{a, t),t). 

Also we define Q{t) = P{^}{t)). Let us note that since i{j and 4> are 2tt periodic, the 
resulting ip and (p are well defined. We do not have problems with the harmonicity of ^/^ or (/> 
at the point which is mapped from minus infinity times i (which belongs to the water region) 
by P since (j) and ip tend to finite limits at minus infinity times i. Also, the periodicity of (p 
and ip causes (p and ip to be continuous (and harmonic) at the interior of P(il^(t)). 

Let us assume that there exists a solution of (LIO) and that we take Unormai = r~^\ such 
that Unormaiioii) , Unormai{o2) < for all < t < T, with T Small enough, thus z{a, t) satisfies 
the arc-chord condition and does not touch the removed branch from P{w). 
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The system (1. 10) in the new coordinates reads 

dni> 



^>„(a,t) 



z(a,t) \Za{a,t)\ 

V = V^'ip in P{n^{t)) 
zt{a, t) = (5^(a, t)u{a, t) + c(a, t)za{a, t) 

$t(a, t) = -Q'^{a, t)\u{a, t)\'^ + c{a, t)u{a, t) ■ Za{ot, t) — gP2^{z{a, t)) 

z{a,0) = 
ic«(a,0) = 4°(a) = $°(a), 



(1.18) 



where u is the hmit of the velocity coming from the fluid region in the tilde domain and 



Q\z{a,t),t) = 



dw 



{p-\z{a,t))) 



, Q\a,t) = 



dP 
dw 



izia,t)) 



We can solve the Neumann problem in the complement of f2(i). Therefore we can represent 
the velocity field v in terms of a vorticity amplitude Q. 
We will see that z and oj satisfy the following equations 



zt{a, t) = Q^{a, t)BR{z, oj){a, t) + c(a, t)za{ce, t). 



(1.19) 



ojt{a,t) = -2dtBR{z,ou){a,t) ■ Za{a,t) - \BR{z,u;)\'^daQ'^ia,t) 

_gj QHa,t) uia,t)^ \ ^2c{a,t)d^BR{~z,u)-U(^,t) 

V 4 \Za(a,t)\^J 

+ da (c(a, t)u{a, t)) - 2gda (iVH^C", *))) ■ (1-20) 

Remark 1.4 Equations (1. 19-1.20) are analogous to (1. 8-1. 9). In fact, if we set Q = 1 in 
(1. 19-1.20) we recover (1. 8-1. 9). 

Our strategy will be the following: we will consider the evolution of the solutions in the 
tilde domain and then see that everything works fine in the original domain. 

We will have to obtain the normal velocity once given the tangential velocity, and vice- 
versa. To do this, we just have to notice that 

^a{oi,t) = u{a,t) ■ Za{a,t) = BR{z,oo) ■ Za{a,t) + ^l^l^l. 

From that, we can invert the equation (see [13]) and get u. Equation (1.6) in the tilde domain 
then tells us v on the boundary dQ,(t). 
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We now note that a solution of the system (1.18) in the tilde domain gives rise to a solution 
of the system (1. 10) in the non-tilde domain, by inverting the map P. In fact, this will be 
the implication used in Theorem I.l (finding a solution in the tilde domain, and therefore in 
the non-tilde). 

Remark 1.5 It is likely that a similar argument works for the other two settings (closed 
contour and asymptotic to horizontal) by choosing an appropriate P{w) that separates the 
singularity. For example, for the closed contour we could consider Pcioi^) = y/z, taking the 
branch so that it separates the singularity, and for the asymptotic to horizontal scenario, it 
is enough to move the interface such that the water region is entirely contained in the lower 

halfplane (and the point —i belongs to the vacuum region) and apply the relation = 

We now state the local existence results that lead to the proof of the existence of a splash 
singularity (Theorem I.l). To avoid the failure of the arc-chord condition, we will prove the 
local existence in the tilde domain. This can be done in two diff^erent settings, namely in the 
space of analytic functions and the Sobolev space H^. 

For the analytic version we define 

Sr = {a + ir], \r]\ < r}, 

_|_ J —IT 

ll/llr = Il/lli2(a5r.) + ll^a/llL2(a5r)' 

we consider the space 

H^{dSr) = {/ analytic in Sr, \\f\\r < oo, / 27r-periodic} 

and we take (zi - a, Z2, $) G {H^{dSr))^ = Xr- 

The first results concerning the Cauchy problem for small data in Sobolev spaces near 
the equilibrium point are due to Craig [18], Nalimov [25] and Yosihara [34]. Beale et al. 
[7] considered the Cauchy problem in the linearized version. For local existence with small 
analytic data see Sulem-Sulem [29]. Our main results regarding local existence in the tilde 
domain are the following theorems: 

Theorem 1.6 (Local existence for analytic initial data in the tilde domain) Letz^{a) 
be a splash curve and let ■ (a) = j^fj- (a) be the initial tangential velocity such that 

(z?(a)-a,zO,(a),$°(«))GXro, 
for some ro > 0, and satisfying: 

^- / «°ormaz(")l^«^°(a)|(^« = 0. 
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Then there exist a finite time T > 0, < r < vq, a time-varying curve z{a, t) and a function 
^{a,t) satisfying: 

1. P~^{z\{a,t)) — a,P''^{z2{oi,t)) are 2iT -periodic, 

2. P~^{z{a,t)) satisfies the arc-chord condition for all t G (0,T], 
and u{a,t) with 

{zi{a,t),Z2{a,t),^{a,t)) G C([0,r],X,) 

which provides a solution of the water wave equations (1.18) with z^(a) = P{z^{a)) and 
u{a,0)-{za)H»,0) = u^{a)-{zX{c^). 

The main tool in the proof is an abstract Cauchy-Kowalewski theorem from [26] and [27]. 
For more details see [11]. 

For the proof of local existence in Sobolev spaces we will take the following c{a,t): 



This choice of c will ensure that \z{a, t)\ depends only on t. We will also define an auxiliary 
function (p{a,t) analogous to the one introduced in [7] (for the linear case) and [4] (nonlinear 
case) which helps us to bound several of the terms that appear: 

^(«,0 = ^'^r^'f ^-ti'^^ - c{a,t)\Uo^,t)\. (1.21) 
2|zQ(a,t)| 

Then, we can prove the following theorem: 

Theorem 1.7 (Local existence for initial data in Sobolev spaces in the tilde domain) 

In the setting of Section LB, let z^{a) be the image of a splash curve by the map P parametrized 
in such a way that \daZ^{a)\ does not depend on a, and such that z^{a), Z2{a) G H'^{T). Let 

(f{a,0) G H^^2(T) be as in (1.21) and let Cd{a,0) G H'^{T). Then there exist a finite time 
T > 0, a time-varying curve z{a,t) G C{[0,T]; H'^) , and functions uj{a,t) G C{[0,T]; H'^) 
and (f G C{[0,T]; H^'^^'j providing a solution of the water wave equations (L19 - L20). 

The proof is based on the adaptation of the local existence proof in [13] to the tilde 
domain. 

Some of the relevant estimates from [13] obviously hold here as well, with essentially 
unchanged proofs. We state such results in Lemmas IV.2 and Lemmas IV.5, . . ., IV.9 below; 
and refer the reader to the relevant sections of [13] for the proofs. 

However, [13] contains several "miracles", i.e., complicated calculations and estimates that 
lead to simple favorable results for no apparent reason. To see that analogous "miracles" occur 
in our present setting, we have to go through the arguments in detail; see Lemmas IV. 10 and 
IV.12, . . ., IV.15 below. 
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We have tried to make it possible to check the correctness of our arguments without 
extreme effort, and without undue repetitions from [13]. 

It would be very interesting to understand a-priori why the "miracles" in this paper and 
in [4], [13] occur. Presumably there is a simple, conceptual explanation, which at present we 
do not know. 

At the end of Section II we will define the notion of a "splat curve" . The curve depicted 
in Figure 2(b) is an example of a splat curve. 

In the statement of Theorem 1.7, we may take zo{a) to be the image of a splat curve 
under P rather than the image of a splash curve. 

The proof of Theorem 1.7 goes through for this case with trivial changes. Consequently, 
we obtain an analogue of Theorem 1.6, with hypothesis 1 replaced by 

Hypothesis I': vf^^^^^^ = Unormal{oc, 0) is negative for all a G /i U /2, where h, h are the 
intervals appearing in the definition of a splat curve in Section II. 

Just as Theorem 1.6 implies the formation of splash singularities for water waves, the 
above analogue of Theorem 1.6 for splat curves implies 

Corollary 1.8 (Splat singularity) There exist solutions of the water wave system that col- 
lapse along an arc in finite time, hut remain otherwise smooth. 

I.D Further Results 

Here we mention some immediate consequences of our results which are relevant: 

1. (Splash and Splat singularities for 3-D water waves) It is possible to extend our re- 
sults to the periodic three dimensional setting by considering scenarios invariant under 
translation in one of the coordinate directions. While preparing the final revisions of 
this manuscript, we noticed that in a very recent arXiv posting [17], Coutand-ShkoUer 
consider additional 3D splash singularities. 

2. (No gravity) The existence of a splash singularity can also be proved in the case where 
the gravity constant g is equal to zero, as long as the Rayleigh- Taylor condition holds. 

II Splash curves: transformation to the tilde domain and back 

In this section wc will rewrite the equations by applying a transformation from the original 
coordinates to new ones which we will denote by tilde. The purpose of this transformation 
is to be able to deal with the failure of the arc-chord condition. 

For initial data we are interested in considering a self- intersecting curve in one point. 
More precisely, we will use as initial data splash curves which are defined this way: 

Definition II.l We say that z{a) = {zi{a), Z2{a)) is a splash curve if 
1. zi(a) — a, Z2{a) are smooth functions and 2Tr -periodic. 
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Figure 4: Tilde domain at times t = (Red - splasli), t = A ■ 10 ^ (Blue - turning) and 
t = 7 ■ 10~3 (Black - graph). 




-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 



X 

Figure 5: Zoom of the splash singularity at times t = (Red - splash), t = A ■ 10~^ (Blue - 
turning) and t = 7 ■ 10~^ (Black - graph). 
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2. z{a) satisfies the arc-chord condition at every point except at ai and a^, with a\ < a.^ 
where z{q.\) = z{a2) and |^;a(ai)|, |2;q,(q;2)| > 0. This means z{ai) = 2(0:2), but if we 
remove either a neighborhood of ai or a neighborhood of 02 in parameter space, then 

the arc-chord condition holds. 

3. The curve z{a) separates the complex plane into two regions; a connected water region 
and a vacuum region (not necessarily connected). The water region contains each point 
x-\-iy for which y is large negative. We choose the parametrization such that the normal 
vector n = (~ ^ ^^"^(q) ^ points to the vacuum region. We regard the interface to be 
part of the water region. 

4. We can choose a branch of the function P on the water region such that the curve 
z{a) = {z\{a) , Z2{a)) = P{z{a)) satisfies: 

(a) zi(q) and Z2{a) are smooth and 2Tr -periodic. 

(b) z is a closed contour. 

(c) z satisfies the arc-chord condition. 

We will choose the branch of the root that produces that 

lim P{x + iy) = -e"^^/'' 
2/— >■— 00 

independently of x. 

5. P{w) is analytic at w and ^{w) if w belongs to the interior of the water region. 
Furthermore, (±7r, 0) and (0, 0) belong to the vacuum region. 

6. z{a) 7^ for Z = 0, 4, where 

^"-^"•"^^ ^'-{t.^t^^ ^'<t.'T^' ''Kt^'v^)' ^'~-{T.'ii}- 

(11.1) 

Prom now on, we will always work with splash curves as initial data unless we say other- 
wise. Condition 6 will be used in the local existence theorems and can be proved to hold for 
short enough time as long as the initial condition satisfies it. It is also immediate to check 
that the previous choice of P transforms any periodic interface into a closed curve. Here are 
two examples of curves which are not splash curves (see Figure 6). 

Now we will show a careful deduction of the equations in the tilde domain. Prom the 
definition of z we have that 



Za{a, t) = VP{z{a, t)) ■ Za{a, t) 



(11.2) 



and 



zt{a, t) = VP{z{a, t)) • Zt{a, t) = VP{z{a, t)) ■ {u{a, t) + c{a, t)za{o(, t)) 
= VP{z{a, t)) ■ u{a, t) + cza{a, t). 



(11.3) 
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Figure 6: Two examples of non splash curves. 

Since cj) = tj) o P and v = Vcf) = V((/) o P), we obtain 

dp 

This implies that 

u{a,t) = VP{z{a,t))^u{a,t). 
Plugging this into (II. 3) we get 

zt{a, t) = VP{z{a, t)) ■ \7P{z{a, t))^ • u{a, t) + cza{a, t). 
From the Cauchy-Riemann equations 



VP{z{a, t)) ■ VP{z{a, t)f = Q'\a, t) ■ /da, Q'\a, t) 



dP{z) 



dz 



In this particular case, this means that 



1 + z(a,t)^ 



5(a, t) = zi{a, t) + iz2{a, t). 



(11.4) 

(II.5) 
(11.6) 

(11.7) 



Az{a,t) 

Recall that ^ is the restriction of (p to the interface, i.e. ^{a,t) = (j){z{a,t),t). Then 



$(a,t) = (l){z{a,t),t) = (t){P~\z{a,t)),t) = 4>{z{a,t),t) = <^{a,t) 



(11.8) 



Thus, $ satisfies 



_ 1 

~dt ~2 
1 



u{a, t)p + c{a, t)u{a, t) ■ Za{a, t) — gz2{a, t) 



\VP{z{a, t)y' • u{a, t)|^ + c(a, t)u{a, t) • Zaia, t) - gP^'^izia, t)), (II.9) 



16 



where the subscript in the gravity term of the last hne denotes the second component. 
Thus the system (1. 10) in the new coordinates reads 

A-i/;{x,y,t) = in P{n'^{t)) 

7 4>a(n.:/) 

2(a,t) \Za.[pL,t)\ 

v = V^i) in P{^'^{t)) 
zt{a, t) = Q^{a, t)u{a, t) + c{a, t)za{a, t) 

$t (a, t) = t)\u{a, t) ^ + c(a, t)u{a, t) ■ z^ia, t) - gP2\z{a, t)) + p*{t) 

z(a,0) = z°(a) 

l'«(a,0) = l>° («) = <(«). (11.10) 
We have seen that v can be represented in the form 



27r J_^\{x- zi{a,t),y- Z2{a,t))\^ 
Taking hmits from the fluid region we obtain 

u{pi,t) = BR{z,Lb) + Za- 

The evolution of u is calculated in the following way. First, let us recall the equations 
zt(a, t) = Q^ia, t)u{a, t) + c(a, t)za{a, t) 

4>t(a,t) = ^Q'^{a,t)\u{a,t)\'^ + c{a,t)u{a,t) ■ Za{a,t) - gP2^{z{a,t)) 

$0(0, t) = u{a^ t) ■ Za{a, t) 
z{a,0) = z°(a) 

0) = $° (a) = <(a). (11.11) 
Substituting the expression for u{a,t) and performing the change c{a,t) = c{a,t) + 
^^^(«'*)p!fe§F obtain 

zt{a, t) = Q^{a, t)BR{z, a))(a, t) + c{a, t)za{a, t) 
^a{a,t) = BR{z,oj){a,t) ■ Za(a,t) + ^w(a,i) 

4>t(Q;, t) = ^(5^(q;, t)\u{a, t)\^ + c(a, t)u{a, t) ■ Za{a, t) - gP^'^{z{a, t)) 
- -Q\a,t)\BR{z,u){a,t)\^ ^ ^ ' ^ ^ ' ^ 



2^ ^ ' ^' V , .V w, 8 |za(a,t)|2 

+ c(Q;,t)Si?(z,a;) • Za{oi,t) + ^c(a, t)a;(a, t) - gP2^{z{a,t)). (11.12) 
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On the one hand, by taking derivatives with respect to t in the second equation follows 

t) = dtBR{z, u;){a, t) ■ Za{a, t) + BR{z, u;){a, t) ■ Zat{a, t) + ^^^^^ 

= dtBR{z,oj) ■ Za{a,t) + \BR{z,u})fdaQ^{a,t) 

+ Q'^{a,t)BR{z,uj) ■ daBR{z,ui) + Ca{a,t)BR{z,oj) ■ Za{a,t) 

+ c{a,t)BR{z,Co) ■ Zaa(.a,t) + (11.13) 
On the other, taking derivatives with respect to a in the third equation in (11.12) yields 

^atia,t) = ^\BR{S,{;j)\'^daQHa,t) + Q\a,t)BR{z,Co) ■ daBR{z,Co) 

~2^"( 4 \~za{at)\y +^c.{a,t)BR{z,Cj)-Za{a,t) 
+ c{a, t)daBR{z, oj) ■ Zaia, t) + c{a, t)BR{z, Cj) ■ Zaa(,<^, t) 
+ Ida {c{a,t)u{a,t)) - da {gP^\z{a,t))) . (11.14) 

Combining both equations, we find that 

Cjt{a,t) = -2dtBR{z,co){a,t) ■ Za{a,t) - \BR{z,oj)fdaQ'^{a,t) 

_^^J^Q2^^^n^ +2c(a,t)a„i?i?(z,a;)-5«(a,i) 

+ da {c{a,t)oj{a,t)) - 2da {gP^\z{a,t))) . (11.15) 

We will proceed in the following way: we will consider the evolution of the solutions in 
the tilde domain and see that everything works fine in the original domain. For example, the 
sign condition on the normal vectors in the non-tilde domain has an equivalent form in the 
tilde domain (i.e. the two normal components have negative sign). 

In the non-tilde domain, this implies that the interface moves away from the branch 
removed from the square root, and therefore the interface touches neither the branch cut nor 
the conflictive points (see Condition 6 in Definition II. 1). Hence P and P~^ will be well 
defined and one-to-one. (See Figure 9). 

Let us note that getting = <^ o P is not a problem since is bounded and harmonic. 
Moreover, as {; = V"'-?/' and 

V = VP'^{v o P) 

and VP has exponential decay at infinity, the velocity v belongs to L'^{Q,'^{t) D [— 7r,7r] x M). 

Remark II.2 ^, u and z have easy transformations to the tilde domain but to has not. 

We would like to discuss what happens to the amplitude of the vorticity w in the non-tilde 
domain as the curve approaches the splash. 

If the vorticity belongs to C{[0,Tspiash],C^{^))-, then the normal velocity should be con- 
tinuous at the splash point and therefore the normal component of the restriction of the 
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velocity to the curve from the water region cannot have the same sign at z{ai) and 2(02) 
(see Theorem I.l). This means that the —norm of the amphtude of the vorticity becomes 
unbounded at the time of the splash. 

We illustrate this phenomenon by plotting 1/ max |a;| (see Figure 7), where the blue curve 
is the calculated cj and the red curve is a potential fitting to the data as numerical instabilities 
don't allow us to compute a; with enough precision when we are in the regime which is close to 
the splash. Time has been reversed so that the splash occurs at time t = and the interface 
separates from itself at t > 0. 




0.5 1 1.5 2 



X 10"* 



Figure 7: Vorticity amplitude in the nontilde domain. The vorticity reaches infinity at a rate 

^ ~ ^ Thefit isgivenby F = 23.72-t°-96'5-1.476-10-^ 



of approximately 



splash' 



-t) 



i^splash ^) 



We also have performed numerical simulations in order to get a blowup rate for the arc- 
chord condition. As in Figure 7, we plot the inverse of the arc-chord constant. The blue 
curve is made by the calculated points and the red curve is the interpolating one. We see a 
very good fitting. Time follows the same convention as before and the numerical evidence 



indicates a blowup of the arc-chord as y 



1 



plash 



The results can be seen in Figure 8. 



,x 10" 




Figure 8: Arc-chord condition in the non-tilde domain. The arc-chord reaches infinity at a 
rate of approximately jTp — ^ . The fit is given by = 11.41 • t + 5.104 • 10~^. 

y-'- splash ^/ 
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We also kept track of the energy conservation. If we consider the following energy (not 
to be confused with the one in Section IV): 

Es(t) = \l \v{x,y,t)fdxdy + l [ giz2{a,t)fdaZi{a,t)da = Ek{t) + Ep{t) (11.16) 

^ Jn'j{t) ^ J—K 

where z{a.,€) = {z\{a,€),Z2{aJ^),u{a^£) = v{z{a,t),t), and = f7^(t) n [— vr,7r] x M 

is a fundamental domain in the water region in a period, then we can see that the energy is 
conserved; this is a check of the accuracy of our numerics. 



dt 



= / v{x,y,t){vt{x,y,t) + v{x,y,t) ■Vv{x,y,t))dxdy 
v{x, y, t){-Vp{x, y, t) - g{0, l))dxdy 



/(*) 



= - / v{x, y, t){V{p{x, y, t) + gy))dxdy 

Ja.){t) 

= - v{x, y, t) ■ itgyds 

Jdin'jit)) 

/TT 
gz2{a,t)u(a,t) ■ daZ^{a,t)da (H-l''') 
-TT 

where we have used the incompressibility of the fluid (V • v = 0) and the continuity of the 
pressure on the interface {p*{t)\QQ2(^f-^ = 0). Next 

dE (t) 1 

— ^— = y gz2{oi,t)dtZ2{oi,t)daZi{a,t)da + - j g(z2{a,t)fdtdaZi{a,t)da 

/TT /"TT 
gz2{a,t)dtZ2{a,t)daZi{a,t)da- / gz2{a,t)daZ2{a,t)dtZi{a,t)da 
-TT J —IT 

/TT 
gz2{a, t)u{a, t) ■ daZ'^{a, t)da. (11.18) 
-TT 

This proves that the energy is constant. Note that: 

/ \v(x,y,t)\'^dxdy = \V4){x,y,t)\'^dxdy 

= -/ (l){x,y,t)A4>{x,y,t)dxdy 
jQ.){t) 

+ I (f){x,y,t)'V(f){x,y,t) -itdxdy 

Jd(n-j(t)) 

^■^"^ / (t){x,y,t)V(l){x,y,t)-ltdxdy (11.19) 
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so the numerical calculation is restricted to the values at the boundary. We observe that 
the energy of our system is conserved, as we have 



max Es{t) — min Es{t) 

Es(t) 38.3936, — ^ 6 • 10"^^ 

minEs{t) 

We now give the proof of Theorem I.l using Theorem 1.7. 

Proof of Theorem I.l: Using the fact that there is local existence to the initial data in the 
tilde domain and applying to the solution obtained there, we can get a curve z(a, t) that 
solves the water wave equation in the non tilde domain. Details on the local existence in the 
tilde domain are shown below. Note that the sign condition (1.12) assumed in Theorem I.l 
guarantees that for positive time t the curve in the nontilde domain will separate (as depicted 
in Figure 9(a)) instead of crossing itself (as depicted in Figure 9(b)). More precisely, we check 
that for small positive time t the curve a ^ z{a,t) = {zi{a,t), Z2{a,t)) = P^^{z{a,t)) G 
M/27rZ X M is a simple closed curve, i.e. that a ^ z{a,t) is one-to-one. Indeed, if not, there 
exist a sequence of positive times ti, and points a^,, such that a^, ^ a^, mod 27rZ, 
but z{aj^,ti,) = z{ai^,tu). Since the initial splash curve a i— t- z{a,0) satisfies the modified 
chord-arc condition described in Condition 2 of Definition II. 1, we may assume without loss of 
generality that — >■ ai and — > 02 (with ai,a2 as in Definition II. 1). The sign condition 
(1.12) therefore guarantees that (for large v), z(a^,ti,) and z{aj^,ty) lie in the image of the 
(open) time-zero water region under the map P. Moreover (for large v), z(a^, in) ^ z{a^^, ti,) 
since z(ai,0) 7^ i(a2,0). 

Since P~^ is one-to-one on the image of the open time-zero water region under P, it follows 
that (for large u) we have z{aj,,ti,) 7^ z{al,ti,) G M/27rZxM, with z{a.,t) = P~^{z{a,t)). This 
contradicts the defining condition z{a^,t,y) = z{aj^, ti,), completing the proof that a i->- z(a, t) 
is a simple closed curve for small positive t. 

The proof of Theorem I.l is complete. □ 

We end this section by defining a "splat curve", as promised in Section I. To do so, 
we simply modify our Definition II. 1 for a splash curve, by replacing Condition 2 in that 
definition by the following 

Condition 2': We arc given two disjoint closed non- degenerate intervals /i,/2 C [0, 27r) 
whose images under a ^ {zi{a), 2:2(0)) G M/27rZ x M coincide. 

The map a >->■ (zi(a), ^2(0)) G M/27rZxM satisfies the chord-arc condition when restricted 
to the complement of any open interval J such that J D /i or J D /2- 

As promised, the curve depicted in Figure 2(b) is a splat curve. Observe that the curve 
in Figure 2(b) cannot be real-analytic. 



Ill Proof of real-analytic short-time existence in tilde domain 

The main goal of this section is to prove Theorem 1.6. In order to accomplish this task we 
will prove local well-posedness for the system (III.l) below. In this section, we will drop the 
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tildes from the notation. The system arises from (11.11) taking c = 0: 

|2 



\dP 
I dw 



U 
2 



l|dP(p-l(^))|2|^|2_^p-l(^^ 



2 I dw 

BR{z,aj) + 



TlZa 



{P-\z{am\ 



2 

_1_ 
16 

log 



2|2a| 

+ BR{z,u;) ■ Za 

zi(a,t)+iz2{a,t) 
i+(zi{a,t)+iz2{a,t))'^ 
i—{zi{a,t)+iz2{a,t)y-' 



(111.1) 



We demand that z°(a) ^ (0,0) to find the function ^{p-'^{z{a,t))) well defined. This 
condition is going to remain true for short time. We also consider z^{a) ^ qi, I = 1, ...,4 in 
(II.l) to get P2^{z{a,t)) well defined. A gain this is going to remain true for short time. 

The main tool in this section is a Cauchy-Kowalewski theorem (see [10, Section 5] for 
more details). We recall the following definitions 

Sr = {a + irj, \r]\ < r}, 



2 

L^{dSr) 



= Y, \f{a±ir)fda, 



Il2(95,.) + ll^a/llL2(a5r)' 

the space 

H^{dSr) = {/ analytic in Sr, \\f\\r <oo,f 27r-periodic} 
and we now take {zi, Z2, $) G {H^{dSr))^ = X^- We have the following theorem: 
Theorem III.l Let z°(a) he a curve satisfying the arc-chord condition 



|z°(a) -z°(a-/3)p 



> 



M2 



which doesn't touch the points qi, I = 0,...,4 in (II.l), and (z , $ ) G X^g for some ro > 0. 
Then, there exist a time T > and < r < ro such that there is a unique solution to the 
system (III.l) in C{[0,T],Xr) with initial conditions z{a,0) = z'^{a),^{a,0) = $°(a), for 
all a G T. 

Equation (III.l) can be extended for complex variables: 

zt{a + i^, t) = F\z{a + i^, t), $(a + i^, t)), $t(a + i^, t) = F\z{a + i^, t), $(a + i^, t)). 
Here 



F\z,<P) 



dP 
dw 



iP-\z)) 



u 
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where we abuse notation by writing 



dP 

dw 



{p-\z{a + ii,t))) 



1 nf^J(zi(a + ii, t)-q{f + {Z2{a + i^, t) - g^) 



J \2l 
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(zi(a + i^,t))2 + (z2(a + ie,t))2 



and 



u{a + i^, t) = BR{z{a + i^, t), u{a + i^, t)) + 



uj{a + i^, t)daz{a + iS^, t) 
2 V {dazi{a + i^, t)y + {daZ2{a + iC, i))2 



with 
1 



PV 



BR{z{a + i^, t),uj{a + i^, t)) = 

{z2ia + i^- 13, t) - Z2{a + i^,t), zi{a + i^,t) - zi{a + i^- P,t)) 
271^ ' Jj{zi{a + itt)-zi{a + iS,-l3,t))^ + {z2ia + i(„t)-Z2ia + iC-/3,t)y 

and oj given imphcitly by 

*a = -+BR{z,Uj) ■ Za- 

We wiU also abuse notation by writing |tt|2 for tt^ + u^) even for complex u = (ui,U2)- The 
operator is given by 



F\z,^) = 



1 



dP 
dw 



{P-\z)) 



\u\^-gP^\z) 



where 



1 4 

P2\z{a + ii, t)) = - Y^i-l)' log[(^i(a + ie, t) - + {z2{a + i^, t) - g^)^]. 

Below we will use a strip of analyticity small enough so that the complex logarithm above is 
continuous. We use the following proposition: 



Proposition III.2 Consider < r < r' and the open set O C X^/ given by: 



O = {{z,<^) e :\\zi\\r',\mr' < R, mf \izi{a + i^) - q{f + {z2ia + i^) - q^fl > R ^, 
1 = 0,.. .,4, inf G{z){a + i^,P) > R-"^} 

/3e[-7r,7r] 

with 



G{z){a + i^,P) = 



{zi{a + iQ - zi{a + i^ - /3))^ + (z2(a + iQ - Z2{a + - ^))^ 

^2 



then the function F = {F^ , F'^) for F : O ^ Xr is a continuous mapping. In addition, there 
is a constant Cr (depending on R only) such that 
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\\F{z,n\r < ^^\\{z.mr' (111.2) 

\\F{z\ _ F{z\ ^^)\\r < ^||(^' -Z\^^- (III.3) 

and 

sup \F\z,^){a + iO-F\z,^){a + i^-P)\<CR\P\ (IIL4) 

ct-\-i^GSr 

for z, z^, G O. 

Proof: First we point out that co is given in term of and z by the implicit equation 

^a = '^+BR{z,Uj)-Za= ^(/ + J)(w). 

It is well known that the operator (/ + J) is invertible on for real functions with mean 
zero (see [13, Section 5] for more details). Writing 

/ , - N / , . X 1 f (z(a±ir) — z(B))-^ ■ Zci(a±ir) ,^ 
io{a±ir) = 2^a{a±ir)-— / , / , .^V 2 MP)dP 

one can find that 

11^11x2(95^) < 2 II $0 11^2(95^) + Cr\\uj\\l2i^qSo) 

(where Cr depends on R) for {z, $) G O. The bound of (/ + J)~^ for real functions yields 
II'^IIl2(95o) ^ 2||(/ + J)"^||i2_j.i2||$Q||i2(95g) < Ci? II $a 11x2(950 • 

Thus 

ll'^llL2(a5^) < C'fl||$a||L2(a5^). 

Analogously, one finds that 

\\d'M\ <CRmr- 
This allows us to assert that oj is at the same level as in terms of derivatives: 



MlLHdSr) + ll^>llL2(a5.) < CRlMr < CRimr'. (111.5) 

Then, inequality (111.2) follows as in [10, Section 6.3]. We will see how to deal with the 
most singular terms. For the first term in the norm, it is easy to find that 

\\F{z,^)\\L2^eSr) < CR\\{z,mr < CR\\{z,^)\\r'. (III.6) 

In order to control the second one, we will show how to deal with F^ as is analogous. 
Here we point out that the functions 

dP 1^ 
— {p-\z{a + t^,t)))\ , P^Hzia + i^,t)) 
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have no loss of derivatives and they are regular as long as {z, $) G O. Therefore, in d^F^ the 
most singular term is given by 



dP 
dw 



iP-\zia + i^,t))) 



as the rest can be estimated in an easier manner (see [13, Section 6.1] as an example with 
more details). Prom the definition it is easy to bound in it remains to 

control d^u in L?'{dSr). To simplify the exposition we ignore the time dependence of the 
functions, we denote 7 = a it ir, 

(^1(7) - ^1(7 - + (^2(7) - ^2(7 - P)? = \z{l) - ^(7 - 

and 

(2:2(7 -P)- ^2(7), ^1(7) - ^1(7 - P)) = - zi-y - P))^. 
Next, we split as follows 

dlu = Ii + l2 + h + h + h + h + l.o.t. 

where l.o.t. denotes lower order terms which can be estimated in an easier manner. We have 

{dlz{^)-dlz{^-P))- 



h 



h 



1 r 

2^^^]-^ \zh)-z{^-P)\i 
= —PV r (^(7)-^(7-/3))' 

7!" J-TT 



Ml - mp, 



h = — PV I 



h 



z{l) - z{j - I3)\^, 
1^(7)5^(7) 



h = ~ 



2 \za{l)\l ' 



2 {\Za{l)\ 



2\2 



■daz{i) ■ d^z{-f) 



and 



For Iq we find 



l d^Ml)daZh) 
Jfi = ; ^^TT, • 



Zah)\t 



1, 



l^6||L2(aS.) ^ oll^a^||L-(S.)( inf G{z){j,f3)\ \\d^ij\\L2(^Qs^) 

/3e[-7r,7r] 



and since (z, $) G O we get 



ll^6||L2(as,) < CR\\dlu:\\L2(aSr) 
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by using Sobolev embedding. A simple application of the Cauchy formula gives 

(J 

which allows us to find 



I 7- II ^ Ilfl2 II 



The bound (III. 5) gives finally 



I r II ^ \\ir.\\ 



In a similar way we obtain 

1 1 -^4 1 1 L2 (95.) + 11-^5 ||L2(a5.) < CR\\dtcZ\\L'^{dSr) < ■:^7^M\r' ■ 

In /a we decompose further: = I^^i + I^^2 where 

^^'^ = j\K{^^P)dloj{^-m, 73,2 = \^^H{dlu){^), 

where H denotes the Hilbert transform and the kernel K is given by 

(^(7)-z(7^/g))^ z^(7) 1 
K7)-z(7-/3)|i |z«(7)|^2tan(/3/2)- 

We can integrate by parts d0{—d^u!{'j — /?)) in Is^i to find 

11-^3,1 ||L2(g5^) < CR\\dlu}\\L2(^gs^) < Cii\\^\\r' 

(see [13, Section 3] for more details). The term 13^2 can be estimated by 

WhMlLHaSr) < CR\\H{diu;)\\L^^eSr) = CRWd^uh^oSr) < ^mr'. 
A similar splitting in /2 = I2 1 + /2 2 with 



71" J-n 



(k«(7)l 

(where A = Hda) gives the kernel L as follows 



^2,2 = -^:^^.«(7)-A(5^)(7), 



L{7,P) ■ (5^.(7) -9„.(7-/?)) = - j^^j;^, I^n^) 
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Heuristically, we regard this operator as no better or no worse than a Hilbert transform of 
d^z. It is easy to prove that 

\\h,i\\L-^{3Sr) < CR\\dlz\\L2^aSr) < CrIMt' 
(see [13, Section 6.1] for more details). The term /2,2 can be bounded as fohows 

WhMlL^idSr) < CR\\A(dlz)\\L2(^Qs^) = Ci? 1 1 5^^11^2(95^) < ^^J^\\z\\r'- 

Analogously, for Ji we find 

\\Ii\\L^{aSr) ^ ^73711^11^'- 

This strategy allows us to deal with d^u and therefore with d^F^. The same applies to d^F'^ 
and we can get finally (III. 2). 
To get (III.3) we write 

where 

J,,{cv) = 2BR{z^,co)-zi 
for £ O and j = 1,2. This implies 

^1-^1 = + BR{z\ o?-u^).zl + BR{z\ u^).zl- BR{z\o^) ■ ^ 

which yields 

{u? - u') = 2(7 + J,2)-i($^ _ $1 ) _ 2(7 + J^^r^BRiz^u') ■ zl - BR{z\u^) ■ 4). 
This helps us to find 

\\u^ - u'W^dSr) + Wdlu;^ - dlco^h^aSr) < C{R)i\\^' - ^% + \\z' - z%). 

We use a decomposition similar to the one used to prove (III.2) which allows us to get finally 
(III.3). Inequality (III.4) follows in an easier manner. □ 

Proof of Theorem III.l: We apply the following result of Nirenberg [26] and Nishida [27]. 

Abstract CaTicliy-Kowalewski Theorem: 
Consider the equation 

= F{u{t)) for |i| < 5 (III.7) 

with initial condition 

u{Q) = «o e Xr^ (III.8) 
For some numbers (7, 7? > 0, assume the following hypothesis: 
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For every pair of numbers r, r' such that < r' < r < rg , -F is a Lipschitz map from 

C 

{u € Xr : ||m — n'^llxr < -R} into Xj.', with Lipschitz constant at most -. Then the 

equation (III. 7) with initial condition (III. 8) has a solution u{t) in C{[—5,5\,Xr) for small 
enough r, (5 > 0. 

The above Abstract Cauchy-Kowalewski Theorem is obviously equivalent to a special case 
of Nishida's Theorem [27], although our notation differs from that of [26]. In place of (III. 7), 

Nirenberg and Nishida treat the more general equation ^ = F{u(t),t). 

The proof of the Abstract Cauchy-Kowalewski Theorem in [26] proceeds by showing that 
the obvious iteration scheme 

u^+\t) = u^+ f F{u^{s))ds 
Jo 

converges in X,. for small enough r (depending on t). 

du 

Our system (III.l) has the form — = F(u) for u = {z, Proposition III. 2 tells us that 

the hypothesis of the Abstract Cauchy-Kowalewski Theorem holds for the system (III.l). In 
particular, for i? > small enough, we obtain the arc-chord condition for every u = {z, <I>) 
such that \\{z, <I>) — (z*^, < for any (arbitrarily small) r > 0. 

Hence, the conclusion of Theorem III.l follows from the Abstract Cauchy-Kowalewski 
Theorem. 

□ 

Proof of Theorem 1. 6: 

Applying Theorem III.l, we obtain a solution of the water wave equation, with the correct 
initial conditions, in the tilde domain. Passing from the tilde domain back to the original 
problem, we obtain a solution of the water wave equations as asserted in Theorem 1.6. 

We have to make sure that, for small positive time, the splash curve evolves as in Figure 
9(a), rather than Figure 9(b). 



(a) Good (b) Bad 

Figure 9: Two different evolutions of the interface. 

This is guaranteed by the hypothesis of Theorem 1.6 regarding the sign of the normal 
component of the initial velocity at the splash point. 

□ 
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IV Proof of short-time existence in Sobolev spaces in the tilde 
domain 



In this section we will show how to obtain a local existence theorem for the water wave 
equations in the tilde domain. The proof is based on energy estimates and uses the fact that 
the Rayleigh-Taylor function is positive. 



IV. A The Rayleigh-Taylor function in the tilde domain 

We begin by recalling the function ip{a,t), which will be studied in detail in Section IV. C 
and in the definition of the Rayleigh-Taylor condition, by the expression 



Q ij{a, t) 
2\za{a,t)\ 



- c\za{a,t)\. 



(IV.l) 



Next we introduce the R-T function: 



a = { BRt{z,u}) + j^BRa{z,cb) ] ■ z^ + 



2 LZa 



? 12 



+ Q 



BR{z,uj) + 



UJ 



2| Zq, I 



if ^ 



(IV.2) 



{VQ)iz)-z^ + 9iVP,-'){z)-~z^. 



This function a coincides with the expression z-^{a, t) ■ Vp{z{a, t),t), where p = po P ^. 
Indeed, it is easy to check that 



dt(f> + 



p-gP-'+p*{t). 



And taking the gradient on the equation (IV. 3) yields 

^t + l (VQ2) |^|2 + q2(~ . ^ _ gyp-\ 



(IV.3) 



(IV.4) 



In addition we know that 



v{z{a,t),t) = BR{z,u){a,t) + -^^^^Za{a,t) 

Z\Zce\Ol, t)\ 



and therefore 
di 



v{z{a,t),t) = dtBR{z,u){a,t) + dt 



oj{a, t) 
2|5a(a,t)|2 



Za{a,t) + 



oj{a, t) 
2|5a(a,t)|2 



(IV.5) 



dtza{a,t). (IV.6) 



On the other hand, by using (IV.4) we have 
d 



dt 



v{z{a,t),t) =dtv{z{a,t),t) + idtz{a,t) ■ V)viz{a,t),t) 

= - o i^Q^) - QHv{z{a, t),t) ■ V)viz{a, t),t) 



Vp{z{a,t),t) - gVP^^{z{a,t)) + {dtz{a,t) ■ V)v{z{a,t),t). (IV.7) 
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Furthermore the equation (11.12) together with (IV.5) gives rise to 

2\zaia,t)\^' 



dtz{a,t) =Q'^v{z{a,t),t) - Q ^^"'j^ Za{a,t) + cza{a,t) 



--Q^v{z{a,t),t) - / f ^^-^^fr?, -c\za{a,t)\ ) Za{a,t). (IV.8) 



Therefore by (IV. 1), we obtain 



dtzia,t) =Q^v{zia,t),t)-^ia,t)pp%. (IV.9) 

\Za{a,t)\ 

By introducing (IV.9) in (IV. 7) we have 

j^v{z{a, t),t) = -^ (VQ2) \v^z{a, t),t)f - Q\v{z{a, t), t) ■ V)v{z{a, t),t) 
- Vpizia, t),t)- gVP^^zia, t)) 

+ Q^{v{z{a,t),t) ■V)v{z{a,t),t) - 0{a,t) .l""^"^'*}.. ■ Vv{z{a,t),t). 

\Za{a,t)\ 

Therefore 

piz{a,t),t) = - \ (VQ2) \v{z{a,t),t)\' " '^(«'*)^^'^^ 

-Vp{z{a,t),t) - gVP^\z{a,t)). (IV.IO) 
Next we take a derivative with respect to a in the equation (IV.5) to get 

dav{z{a,t),t) = daBR{z,u)){a,t) + da (^ 2\^^{a^t)\'^ ^ + (^ 21^^(0 ^t)p ) 

(IV.ll) 

Multiplying equation (IV.IO) by Za{(X,t) and using (IV.ll) we learn 
j^v{zia,t),t)^ ■ z^{a,t) = - QVQ ■ z^(a,t)|i)(z(a,i),i)l' 

- ^^^9„i?i?(z,cD)(a,t) • 4(a,t) 

^ ^~za.{a,t).zt{a,t) 



Za 

{a,t)\ \2|^Q,(a, i)p, 
- Vp{2{a, t),t) ■ z^{a, t) ~ gVP^\z{a, t)) ■ z^{a, t). (IV.12) 



On the other hand, by multiplying (IV.6) by z^(a,i) we have 



v{z{a,t),t) ■ Za{a,t) =dtBR{z,ui) ■ z„(a,i) + — — — —7^dtZaia,t) ■ z^ia^t). 
at J 2\Za{oi,t)\^ 



(IV.13) 
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Prom (IV. 12) and (IV. 13) we find 
dtBR{z,co)-z^{a,t) + 



-QVQ-5„(a,t)|v(z(a,t),i)| 



2|za(a,i)P 

2 <^(a,*) 



dtZa{a,t) ■ z^ia,t) 



\zaia,t)\ 



daBR{z,Ld){a,t) ■ z^{a,t) 



(p{a,t) f Cj{a,t) 



a 



\za{a,t)\ \2\zaia,t)\'^ ^ 
-Vp{z{a,t),t) ■ z^{a,t) - gVP^\z{a,t)) ■ z^{a,t). 



(IV. 14) 



Finally, rearranging the terms in (IV. 14) yields 

- Vpiz{a, t),t)- ia,t)= (dtBR{z, uj){a,t) + 



\Zai0!,t)\ 



daBR{z,u}){a,t) ■z^{a,t) 



+ 



uj{a, t) 
2|Sa(a,0P 



dtZa{a,t) + jP^^^Zaa] ' z^ {a,t) + gV P^^ ■ z^{a,t) 



\zaia,t)\ ' 



+ Q 



BR{z,uj){a,t) + 



oj{a, t) 



I Za (ck) ^) 



VQ-z^{a,t)) 



2|z«(a,t)|2' 

and then, comparing with (IV.2), we obtain the desired result 

-Vpizia, t), t) ■ z^{a, t) = (7(a, t). 

Note that for the tilde domain, the Rayleigh- Taylor condition is the same as in the first 
domain, i.e: 

Vp{a,t) ■ z^{a,t) = Vp{a,t) ■ z^{a,t) 

where p = po and 

Zaia, t) = VP{z{a, t)) ■ Za{a, t) =^ z^{a, t) = {-J\7P{z{a, t))J) ■ z^{a, t) 

-1 



where J is the rotation matrix 
this implies that 



1 



. Together with the Cauchy-Riemann equations 



Moreover 



{-JVP{zia,t))J) = VP(z(a,t)). 
Vp{a,t) = VP{z{a,t)fVp{a,t). 



Hence 



{Vp{a,t),z^{a,t)) = {VP{z{a,t)fVpia,t),{VP{z{a,t)))-'z^{a,t)) (IV.15) 
= {Vp{a,t),z^{a,t)). (IV.16) 
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By taking the divergence on the Euler equation (1. 1-1.2) and because the flow is irrotational 
in the interior of the regions Cl^ {t) follows 

-Ap= \Vvf > 

which, together with the fact that the pressure is zero on the interface and p{x, y, t) + gy = 
0(1) when y tends to — oo,then follows by Hopf's lemma in 0?{t) that 

a{a,t) = -\z-^{a,t)\dnp{z{a,t),t) > 0, 

except in the case v = 0. This argument was suggested by Hou and Caflisch (see [31]), 
although the proof of the positivity of the Rayleigh- Taylor condition in the nontilde domain 
for all time was first introduced by Wu in [30] . 

The above proof shows that a > provided our domain (l{t) arises by applying the map 
P to a domain Q{t) with smooth boundary. Here, dQ{t) may be a splash curve, but we 
cannot allow boundaries dO,{t) whose inverse images under P look like figure 9(b). 

Nevertheless, since cr > for the image of P applied to a splash curve, we know that 
(T > at time t = in the context of Theorem 1.7. Our estimates below will guarantee that 
the condition a > persists for a short time. Thus, in proving Theorem 1.7, we may use the 
positivity of a. 



IV. B Definition of c in the tilde domain 



Prom now on, we will drop the tildes from the notation for simplicity. We will choose the 
following tangential term: 



Here and in (n.l2) we find 

i + {zi{a,t) +iz2{a,t))'^ 



(IV.17) 



and 



P^^ = P^\z{a,t)) = log 

Q = Q{z{a,t)) = 



i - {z]_{a,t) + iz2{a,t)Y 

1 + {zi{a,t) +iz2{a,t)Y 



zi{a, t) + iz2{a, t) 

These functions are regular as long as z{a,t) ^ qK We deal with initial data which satisfy 
the above condition and we will show that it's going to remain true for short time. In order 
to measure it we define 

m{q'')(t) = min \ z{a, t) — g'| 

for I = 0,...,4. 

We also point out that, because of our choice of c{a,t), solutions of (1.19 - 1.20) satisfy 
that 

|za(a, i)p = A{t) for any a G T 
as in [14, Equations (2.2 - 2.5)]. 
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IV. C Time evolution of the function in the tilde domain 

Recall that we have defined an auxiliary function t) adapted to the tilde domain, which 
helps us to bound several of the terms that appear: 

yK^) = 01 +\\ c{a,t)\za,{a,t)\. (IV.18) 

We will show how to find the evolution equation for ipt. We have 



and therefore 



that yields 



Q2 A\za? «^ 



The equation for ujt reads: 

ujt = -2BRt • Za - 2QQa\BR\^ + 2cBR^ ■ zg -dg +dg -25a (5^2"') • 



(la) 



For the quantity (1) = (la) + (16) we write 



(i&) 

(IV.19) 



(1) = (la) + (16) = 2cBRg -Zg + dg (^^) = 2c{BRg . zg + - ^M!^) 

|2 /-TT ^„ f/^2T}n\ |2/ 



and then (IV.19) becomes 

/ 2 

a;t = - 2BRt ■ Zg - 2QQg\BR\'^ - dg ' 

|2 f-TT ^„ DD . -V 0^2u |2 



^ ' (IV.20) 



I Ckal r (n^T,T.^ ,0 ^CQgBR ■ Zg 2c'\Zg\'Qg . 
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Furthermore 

UJ 



= QQt 
= QQt 



Za ■ Zat + 1 - Ot{c\Za\ 

^\Za\ 



2\Zaf 

Q'^oj 1 



+ 



Q 



Za 

Za\ 2\Za\ 27r 
2 



{Q^BR)i^ 



2\Za\ 



-2BRfZa-2QQ^\BR\'-d, 



+ 



CiZn 



{Q^BR)p.j-^dp- 



Q 

AcQaBR ■ z, 
Q 



- dt{c\Za\ 



We should remark that we have used that 

1 

Za- Zat= ^ J {Q^BR)p ■ zpdl3. 

For simphcity, we denote 
Computing 

= QQt] — I -w, — ^B{t) - I — -BRt ■ Za : — 7-\BR\ - - — -da I 

\ Za \ ^ \ Za \ \ Za \ \ Za \ ^ \ Za \ \ ^ 

^ V ' 

(2a) 

+c\za\B{t) -^f^BR . Za - - {gP,-') - dt{c\z, 

V \Za\ \Za\ 

(26) 

We can write 

(2) =(2a) + (26) = -B{t)<p, 

and it yields 



+ QQt A - 2c5i? • ^QQa - ^C^\Za\ - ^Qa\BRf - dt{c\Za\). 

I Za I I Za I | Za \ 

We will use the equation above to perform energy estimates. 



(IV.22) 
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IV.D Definition and a priori estimates of the energy in the tilde domain 

Let us consider for > 4 the following definition of energy E{t): 



Eit) = l + \\z\\l,.,it)+ r 9(^\dk,\2da+\\Tiz)\\looit) 



+ Ml,.,{t) + Ml,_,^ it) + ^^^^^^ + E ^^^^y 

where 

\z{a) - z{a - p)\ 

and m{Q'^a) = mma£T{Q'^{z{a,t))a{a,t)}. In the next section we shall show a proof of the 
following lemma. 

Lemma IV. 1 Let z{a,t) and Ld{a,t) be a solution of (1. 19-1.20). Then, the following a 
priori estimate holds: 

j^E{t) < CEP{t), (IV.24) 

for k > A and C and p constants depending only on k. 

The following subsections are devoted to proving Lemma IV. 1 by showing the regularity 
of the different elements involved in the problem: the Birkhoff-Rott integral, zt{a, t), oot{a, t), 
uj{a,t); BRt{a,t), the R-T function a{a,t) and its time derivative at{a,t). 

IV.D.l Estimates for BR 

In this section we show that the Birkhoff-Rott integral is as regular as daZ. 
Lemma IV.2 The following estimate holds 

\\BRiz,u)\\H, < c{\\:f{z)\\1^ + \\z\\l,+, + \Ml,y, (iv.25) 

for k >2, where C and j are constants independent of z and co. 

Remark IV.3 Using this estimate for k = 2 we find easily that 

\\daBR{z,io)\\L^ < C(||J^(z)||ioo + \\z\\ls + Mll^y, (IV.26) 

which shall be used throughout the paper, where C and j are universal constants. 

Proof: The proof can be done as in [13, Section 6.1] since the definition for the Birkhoff-Rott 
operator is independent of the domain. □ 



(IV.23) 
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IV. D. 2 Estimates for zt 

In this section we show that zt is as regular as daZ. 
Lemma IV.4 The following estimate holds 

\\zt\U < c (^\Hz)\\U + MW + Ml. + iz ' (^^-27) 

for k>2, where C and j are constants that depend only on k. 

Proof: It follows from [13, Section 6.2]. The only additional thing we need to control is an 
L°° norm of Q^, which we can easily bound by the m(g') terms which control the distance 
from the curve to the points, more precisely, the one that controls the distance from the 
origin. □ 

IV. D. 3 Estimates for cjt 

This section is devoted to showing that cot is as regular as daOJ. 
Lemma IV.5 The following estimate holds 

MWh'^ < C I ||^(2;)||ioo + ||z||^fe+2 + ||w||^fe+i + 



,0 m{q^){t)^ 



for k>l, where C and j are constants that depend only on k. 



Proof: We use formula (IV. 20) and proceed as in [13, Section 6.3]. Note that in [13] an 
exponential growth appears in the bound of the estimates for the nonlocal operator acting 
on ojt (see equation (IV.20)). However, in a recent paper [15] the authors get a polynomial 
growth for the operator in both 2 and 3 dimensions. Note that even the exponential growth 
is still good enough to prove Theorem 1.7. 

□ 

IV. D. 4 Estimates for u 

In this section we show that the amplitude of the vorticity uj lies at the same level as d^z. We 
shall consider z G H^{T), ip G //^~2 (T) and u G ii'*^~^(T) as part of the energy estimates. 
The inequality below yields oj G H^~^i^). 

Lemma IV.6 The following estimate holds 

< c (^\:f{z)\\1^ + \\z\\l, + Ml,., + Ml>^-^ + E ;;^(^) ' ^^^-^^^ 

for k > 3, where C and j are constants that depend only on k. 
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Proof: We can apply the same techniques as in [13, Section 6.4] since the most singular 
terms are treated there and the other terms are harmless and can be easily estimated. The 
impact of Q is now taken into account by the m{q^) terms (which now cover all of the points 
q^,...,q% □ 

IV.D.5 Estimates for BRf 

Here we prove that the time derivative of the Birkhoff-Rott integral is at the same level as 
dlz. 

Lemma IV. 7 The following estimate holds 

\\BRt\\H^ < C I ||^(2;)||icx, + ||^;||^fc+2 + ||a;||^fc+i + ||(^||^fc+i 



for k >2, where C and j are constants that depend only on k. 



Proof: We proceed as in [13, Section 6.5], where BRt appears in the formula (IV.2). We use 
(IV. 27) and (IV. 28) to bound zt and oot in BRt respectively. □ 

IV. D. 6 Estimates for the Rayleigh- Taylor function a 

Here we prove that the Rayleigh- Taylor function is at the same level as d'^z. 
Lemma IV.8 The following estimate holds 

for k>2, where C and j are constants that depend only on k. 

Proof: We proceed as in [13, Section 6.5] using formula (IV.2). There is a new term in the 

2 

definition of a, namely Q BR{z,uj) + j^^Za {VQ){z) ■ z^, but this term is less singular 
than BRt{z,u) ■ z^. Hence, the new term causes no trouble. □ 

IV. D. 7 Estimates for at 

In this section we obtain an upper bound for the L°° norm of at that will be used in the 
energy inequalities and in the treatment of the Rayleigh- Taylor condition. 

Lemma IV.9 The following estimate holds 

Wath'- < c [wmwl^ + 11^11^4 + Ijcll^a + \Ml^ + J2 ' (1^-32) 



m{q^){t)^ 



where C and j are universal constants. 



Proof: Again, as in the previous subsection, the new term is less singular than the terms 
treated in [13, Section 6.6]. Hence we deal with them with no problem. □ 
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IV. D. 8 Energy estimates on the curve 

In this section we give the proof of the following lemma when, again, k = A. The case k > A 
is left to the reader. Regarding ||5a2;|||2 let us remark that we have 

Lemma IV. 10 Let z{a,t) and u{a,t) be a solution of (1. 19-1.20). Then, the following a 
priori estimate holds: 



f^^\dU'da+\\:F{z)\\l^)<S{t) + CEP{t), (IV.34) 



-IT 

for 



/7f _L 
2Q^a f " 3" A{d^-^^)da, (IV.35) 

and k > A, where C and p are constants that depend only on k. 

(The term S (t) is uncontrolled but it will appear in the equation of the evolution of tp with 
the opposite sign.) 

Proof: Using (IV.27) and (IV. 33) one gets easily 

d 



^^,,^11^3 < C I {\z{a)\\zt{a)\ + \d-^z{a)\\d-^zt{a)\)da 



< CEP{t). 

We obtain 



^^mz)\\l^<CEPit) 



T+ I I — \2\^ocA da= / I I — \d^z\ da+ / - — r^d^z ■ d^ztda 



in a similar manner as in [13, Section 7.2]. It remains to deal with the quantity 

=h+h. 

The bounds (IV.27), (IV.33) and (IV.32) give us 

h < CEP{t). 

Next for I2 we write 

2gV. 



-IT 1 
Jl + J2. 



^Xz-dtiQ''BR)da+ r ^^d^z ■ di{cdaz)da 

J —TT \^Ol\ 
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The most singular terms in Ji are given by Ki, K2, K3 and K4: 



1 r r Q'a ^ {djz - djz') 

K\ = -PV / / I — 12 C'a^ • 1 7j2 ^ dpda, 



= --PV r r ^dtz ■ - z') . {d^z - dtz'Wdpda, 

TT ./_7r J— TT \Za\ \Z Z \ 



^ J —TT J —n \Za\^ \Z Z 

and 



Ki = - r ^d^z ■ BRVQ{z) ■ d^zda, 

^ J —TT I Za I 



where the prime denotes a function in the variable a — /3, i.e. f = f{a — P). 
Then we write: 



That is, we have performed a manipulation in Ki, allowing us to show that Li, its most 
singular term, vanishes: 



= 0. 

The term L2 involves a S.I.O. (Singular Integral Operator) acting on d^z{a) thanks to the 
minus sign between the two terms Q'^au. One can show that 

L2 < C||J^(z)||i.o 11^11^3 ||a;||cM|k|bM||Q^|bM||aM|i. < CEP{t). 

Inside K2 we find that {z — z') ■ {d^z — d^z') can be written as follows: 

{z - z') ■ {dU - diz') = {z-z'- zJ3) ■ {dU - dW) 

- - 4) • diz' (IV.36) 
+ !3{zo^-diz-z'^-diz'), 
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then using that 

za ■ diz = -3dlz ■ dlz, (IV.37) 

we can spht K2 as a sum of S.I.O.s operating on d^z{a), plus a kernel of the form ^^p^ 
acting on d^z ■ d^z with r) e allowing us to obtain again the estimate 

K2 < CEP{t). 

Note that below we will also use a variant of (IV.37), namely 

za ■ {dtz - diz') = (4 - za) ■ diz' - ^{dlz ■ dlz - dlz' ■ dlz'). 
The term is a sum of 



(IV.38) 



L3 = 



12 " 



Z ■ 



\z-z'\'^ |2a|22tan(/3/2)_ 



d^co'd^da, 



plus the following term: 



L4 



/TT 
-TT 



O^z ' z^ 



4 



H{d^uj)da. 



We can integrate by parts on L3 with respect to /3 since d^u' = —di3{d^Lo'). This calculation 
gives a S.I.O. acting on d^oj which can be estimated as before. 
Next in L4 we write 



L4 



and decompose further 

L4 = r 2qv^3 

+ r 2Q^a^^^^^A{d^ip)da 

J —TT 



^ I ^OL I ^ 



da 



J — TT I Zo. I 

= M_i + 5 + Mi, 



for S(t) given by (IV.35). In M_i we find a commutator that allows us to obtain 

M_i < CEP{t). 

Using (IV. 17) for M\ we have 



Ml = -2 / if(gV 



2 d^z -^0x04 



)di{c\za\))da = Ni+N2 + N3 + N^, 
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where 

Ni = 2 H{Q^a^f^)Q^BR^ ■ ^da, 
N2 = 2 r H{Q^a^^^)Q''diBR • ^da, 

J—n \^a\ \^a\ 

J —IT l^ttl \z<y\ 

and A^4 is given by the rest of the terms which can be controhed easily by the estimates from 
Section IV.D.l for the Birkhoff-Rott integral. 

Regarding Ni a straightforward calculation gives 

Ni < CEP it), 

and analogously for 

iVa < CEP it). 

Again, in N2 we consider the most singular terms given by 

02 = - r H{Q'a^^^)^^ . f PF r ^^^(. - z') . {d^z - dU'Wdadp, 

J —TT \^ct\ \Za\ J —TV \^ ^ \ 

03 = 2 r H{Q'a^^^)Q^^ ■ BR{z,dtu)da. 

J —IT I ■^CK I \^Ol\ 

Using the decomposition (IV.36) we can easily estimate O2 as in our discussion of K2. 
In O3 we find 

■ BR{z, dtc.) = I^-J[ ^'-'^'-J^f^^ dtu^'dp. 

Above we can integrate by parts as in our discussion of L3. We find that 

O3 < CEPit). 

Next we split Oi into a S.I.O. acting on (d^z)-^, which can be estimated as before, plus 

the term 



r H{Q^a^f^)Q''ujj% ■ A((a^)^)da. 

J —IT I ^oi I I Za I 



Then the following estimate for the commutator 

||q2^^ . A((9^)^) - A(Q2^^ • idU)^)\\L2 < CEP{t) 
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yields 
where 

Using that 

We can write 



Pi < CEP{t) + R 
R = - r Q^a^^^da{Q^u;j^{d^,z)^)da. 

J — TT \^Ct\ \^Oc\ 

/TT 
f (a) dag {(x) da. 
-TT 



J —n \^ci\ pal J —n \^oi\ 

and a straightforward integration by parts let us control R. This calculation allows us to get 

Pi < CEP{t). 

We can easily show that 

Ki < CEP{t) 

because we can bound Q^aBRVQ in L°°. So finally we have controlled Ji in the following 



manner: 



Ji < CEP{t) + S. 

To finish the proof let us observe that the term J2 can be estimated integrating by parts, 
using the identity d^z ■ daZ = —^d^z ■ d^z to treat its most singular component. We have 
obtained ^ 

/ ^d^z ■ dazd^cda = 3 I y^daiQ^ad^z ■ dlz)dlcda 
and this yields the desired control. □ 

IV. D. 9 Energy estimates for ui 

In this section we show the following result. 

Lemma IV.ll Let z{a,t) and uj{a,t) be a solution of (1. 19-1. 20). Then, the following a 
priori estimate holds: 

j^\M]ju-.{t)<CEP{t), (IV.39) 
for k>A, where C and p are constants that depend only on k. 
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Proof: We will discuss the case k = 4, leaving the other cases to the reader. Formula (IV.28) 
shows easily that 



dt 



Mh{t) < (^\:F{z)\\l^{t) + \\z\\j,,{t) + Mls{t) + MUt) + E 



which together with (IV. 29) yields 

dt 



^ "^\\H2{t) < CEP{t). 



□ 



IV.D.IO Finding the Rayleigh- Taylor function in the equation for da^t- 
In this section we get the R-T function in the evolution equation for da<pt- 

Lemma IV. 12 Let z{a,t) and Lo{a,t) be a solution of (1. 19-1. 20). Then, the following 
identity holds: 

^o.t = NICE - ^^'aa - QV^f^ (IV.40) 
I I I I 

where NICE satisfies 

r A{d^-^cp)d^-^{NICE)da < CEP{t) (IV.41) 

J —IT 

and 

\\NICE\\Hk-2 < CEP{t) (IV.42) 
for k > A, where C and p are constants that depend only on k. 

Proof: We will give the proof for k = 4. Prom now on, when we show that a term / satisfies 
r K{dl^)dlfda < CEP{t) and < CEP{t) 

J —IT 

we say that this term is "NICE" . Then, / becomes part of NICE and by abuse of notation 
we denote / by NICE. Notice that, whenever we can estimate the norm of A^^'^d^f by 
CEP{t), then / is NICE. 

We use (IV.22) to compute 

- ^ - (1^ (I) J 

{3a) 



\ pal / n \ ^ /« \ 



a \ /a \ I "I / a 



-{c\Za\)at ■ 

' 

(36) 
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Expanding (3) = (3a) + (36): 



(3) =(3a) + (36) = - Q'BRt • ^ - {c\za\) 



jat 



\Za\ \pa|/a V / t 

2dD I ^'^ \ l\ \TD(4-\\ I //^2 D D\ I ^'^ \ I O/zO/n DE>\ ^'^ 



= -Q'BRf [j:^] -{\za\B{t)\ + iQ'BR)a- [y^] +2{QQtBR)a- 
We use that 



Za \ Z(yc^ ' Z^ Z^ 1 Zqi \ Z(yt ' Z^ Z^ 



^ot \ J n I I ^ I I V I \ J f I I ^ I I 



to find 



(4) 



^\Za\ Vpal^; / \Za\ 



(5) (6) 



+ {Q'BR)a . zt^^^ + 2{QQtBR\ • ^ - Q^ ^" \ \ '"^ +[QQt 



Za I I I 



(13) (7) (8) (9) 

-{2cBR.^<QQ^ -[^c^\zS^ -{£-\,BR\''q}j . (IV.43) 
V 2- V ^ 3 V ^ 2- 

(10) (11) (12) 

The term {\za\B{t))t depends only on t so it is going to be part of NICE. 

(4) = -B{t)ipa is NICE (at the level of ifa). 

(5) = ^1 p = -] 1 - ] 

Z I Za I I Za I I Za \ 

The first term is at the level of (pa so it is NICE. The second one is the transport term which 
appears in (IV. 54). 

^ ' HkalQ^y \Za\Q^ \Za\Q ^ Q \\Za\ 

Above we find the first term at the level of Za so it is NICE. The second term is at the level 
of (ypQ so it is NICE. Wc write the last one as 



The first term is at the level of Za so it is NICE. For the second term we have used that 



44 



Finally: 



(6) = NICE +^VQ-4^f^. 
(7) = 2{QQtBR)a ■ ^ = 2Q^QtBR ■^+2Q ( ^) BR-Za + 2QQtBRa ■ 



I ^oi I I ^oi I \ I \ / (X I '^Q I 



The first term is at the level of Za^zt^ BR ~ Za so it is NICE. We use that 
Using equation (IV. 21) 



Za ■ Zat ^ ^^^^ 



I ^oi 



and 



we find that 



^=zt.{ V^Q{z) . + VQ(z) • z^'-^ 



\z, 



+ VQ{z)-^B{t). 



That yields 



(7) = 2{QQtBR)a^ = NICE + 2QBR -ZaZf (V^Qiz) ■ ^ 

Zfy \ Zn 



NICE (at the level of Zc,,zt,BR) 



+ 2QBR ■ ZaVQiz) ■ z^ ^"*" 3" + 2QBR ■ Zo,VQ{z) ■ ^S(t) 



+ 2QQtBRa 

\Za 



NICE (at the level of Za,zt,BK) 

Zqi 



(IV.44) 



Finally: 

(7) = 2{QQtBR)a^^ = NICE + 2QBR ■ z«VQ(z) • 4 ^°*",^" + 2QQtBRa ■ 
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\ pal /a \ l^alJa 

= - 2QVQ ■ ZaVgP^\z) ■ ^ - Q^z^ ■ f V V2"'(-^) • 

I ^o: I \ I I / 

NICE (at the level of Za) NICE (at the level of Za) 



-Q^VgP^\z)-z, 

which means 



\z |3 



Next 



(8) = - fQ^i9^\z)\,\ ^ ^^^^ _ Q2^gp-l^,^ . ^± 



(9) = { gg* A1 = QaQ* A +Q?^,oj ^qqJ 



I Zn I ^ I ■^q; I / a 



NICE (at the level of Za,zt) 

We use (IV.44) to deal with We find that 

(9) = [qQi^^ = NICE + QuVQ{z) ■ + QQt f ^ 

For the next term 

(10) = - (2cBR ■ ^QQa) = -2cBR ■ ^,QI - (2cBR ■ ^ J QQo, 

V Fa I Jo. ^ pal V pal/ a 

NICE as before 

- 2cBR ■ ZaQVQiz) ■ -ga ^T 'ia" -^cBR ■ ^Qz^ ■ (V^QC^)) • • 



NICE as before 

Therefore 



(10) = - 2cSii: • p^QQa = NICE - 2cBR ■ ] QQ^ 

\ pal J a \ pal 



- 2cBR ■ ZaQVQ ■ Z, 



U 13 



46 



Next 

The fact that the last two terms are NICE, allows us to find that 



(11) = - d^^'l^-l)^ = NICE - {c'\zg\)^ ^ - ^VQ(.) • 4^ 
Finally: 

m = - ( ^\BR\^Qg] =-^-^^$^-^i\BR\XQg 



Zg / „ Zg Zg I 



which implies that 



NICE 

-^\BR\'zg ■ (V'Qiz) ■ Zg) -Q'\BR\'VQiz) ■ 

Zg Zg 



NICE 



(12) = - (^\BR\'Qg) = NICE - ^{\BR\XQg - Q-'\BR\^VQ{z) ■ z^^-^f^ 

\ I ^Ol I J g I ^Ol I I '^Ol I 

We gather all the formulas from (4) to (12), keeping term (13) unchanged. They yield: 

09 

(fgt =NICE - -M)gg 

\Zg\ 

2 J_ J_ _L 

+ -y^^Qiz) ■ z^ — — Q BRt ■ z^ — — Q VgP2 {z) ■ z^ 



. ~g I |3 ^ ^^'^i -a i 13 ^ - i)^ 2 -a 1 13 

V '"^ V '"^ V 

(I60) (15a) (156) 



+ QuVQiz) ■ 4 + QQt ( ^ ) + 2gBi? • ZgVQ{z) ■ z^ 



(18a) (14a) (186) 

+ QQt2BRg ■ ^ - (2cBR ■ ^) QQg -2cBR ■ ZgQVQ{z) ■ z^ 

\Zg\ V 

^ V V 

(146) (17a) (166) 



\z |3 
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(176) (16c) (17c) 



(led) 

We compute 

(14) = (14a) + (146) = QQt (^) + QQt2BR^ ■ ^ 

= - 2|(Q^)„^ - ^'^(QXBR . ^ + 2|(|.„|BW). 

The last formula allows us to conclude that (14)=NICE. 
We reorganize gathering 

(15) = (15a) + (156), 

(16) = (16a) + (166) + (16c) + (16d), 
(17) = (17a) + (176) + (17c) 

and 

(18) = (18a) + (186) 

as follows: 

I 



ifat =NICE - y^ipaa -Q^BRt ■ + VgP^\z) ■ 



I I I I ^ 

V ' 

(15) 

" V ' 

(16) 

+ iQ^BR)a ■ z^'f^ + {Qu + 2QBR • Za)VQiz) • z^f^ 

\Za\ Za 



(18) 



(17) 
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We add and subtract terms in order to find the R-T condition. We recall here that 

('■P \ ± f V \ 
BRt + I .BRa j ■ + 12 I Zat + ] r -Zaa 1 



■cto. I 



+ Q 

Then, we find 



iPat =NICE - I -ifaa 



BR + —. r^Za 

^\Zar 



VQ{z)-z^+VgP^\z)-zi 



((BRt + j^BRa) ■ + (z^t + ^^aa) • z^ + VgP^\z) ■ z^\ ^ 

V V \Za\ J '^\Za\ \ \Za\ J J \Za 



3 



(19) 



a 



+ (Qui + 2QBR ■ 0„) VQ(i) ^ zi '"^^ 



Line (19) can be written as 



(19) =(Q BR)a ■ Z^ + Q BRa ■ Z^- r— p^— + Z^t + :Zaa ' Z^ 



3 ' """^ ° Ir I Ir |3 ' OU 12 V^"': ' U -a | 13 

(Q BR)a ■ Z^ + (Q BR)a ■ Z^j— + Zat ■ Z^ + J—^Zaa ■ Z^ 

\Za\ \Za\ \Za\ ^\Za\ \ \^ct\ 

^ Zaa ' Z(y 



2QQaBR-z^ 



I Za I I Za I 



Zqi \ ^ 1 3 I Zat ' Za ~\~ , , Zaa ' Za \ ^ r\\ 1 2 I 13 1 ' ^O- ~^ I „ I ^cioi ' Za ) Z^ 



Za Za I 



(zat ' I ' (^{Q BR)a ■ Za + ^2^°"^ ' ^ 



2QQaBR.z^r^'-^f^ 

I Za I I •^a I 
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We expand Zat to find 



(19) 



iQ''BR)a ■ 4 + 



2\zc 



1 2 ^otct ' ^oc 



I I I I 



We denote 

We claim that 

where 

That means 
Thus 

We write 

D{a) 



D{a) = iQ^BR)o, ■ 4 + 



D{a) = AN(a) + \za\H{daip) 
\\AN\\hs < CEP{t). 
{D{a)f = NICE. 



(19) = NICE - 2QQaBR ■ z^ 



\^a\ I I 



2QQ„BR-z^ +Q^^PV 



2- 



Z — 



part of AN, at the level of Za ^ 



-Q'-PV I ^\ ''^'"% z - z') . (.„ - z'^)Jdp 

TT ./ \Z — Z \^ 



part of AN, we use (IV.38) 

.1 \, J . 



part of AN, we use (IV.36) and (IV.37) 



+ Q'^BR{z,uja) ■z^ + 

^ V ' 

AN + ^H{oJc) 



Zrsn ' Z^ . 



Therefore 



D{a) = AN + \za\Q^H 



AN + \zc\H 



UJ 



2bn 



'2\^a\JJ 2| 



+ 



1 2 ■^ao; ■ 



AN+|z«|//(5„<^) + //((c| 

Za\ )a) ~l~ 



2ka 



2 -^ao: ■ 



AN + \za\H{^a) - H {{Q^BR)a ■ Za) + 



2\zo 



2 Zaa ' . 



(IV.45) 

(IV.46) 
(IV.47) 
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We have 



AN 

-Q'-PV [ ; {z - z') . {z^ - z'^Wdp 

TT J \z — zy 



AN, we use that {^z—z')^-Za=(z—z'—fiza)^- 



AN, we use that (z—z')^-Za=(z—z'—^Za)^-Za 
For the second term on the right one finds 



\z-z'\^ 2tt J \z-z'\^ ^ 

2TT J \Z-Z'Y ^-K J \Z-Z'Y 

- ^PV J ^"""."^^/u" ''" (^ - • (^a^ - dlz'P'dp + Lo.t., 



where in Let. we gather the terms of lower order. Then, all the terms above can be estimated 
in but the first one on the right. That is equal to 

plus a commutator which can be estimated in L^. This means that 

{Q''BR)a • z« = AN + iif (q^^f^u 
Taking Hilbert transforms: 



-H {{Q'BR)^ . z^) = AN - iq^'^Jfl^uX = AN + lg2^ia_i«^_ 

^ \ I Za I / -^1 Zn I 



Using that z^^ ■ Za = —Zaa ■ Za we complete the proof of (IV. 46). Thus (19) yields 
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if at =NICE - -Apaa 

I I 



-Q3(|B«P + ^+2.^-|,)vO(.)..i 
+ (Qu + 1QBR ■ Zc)VQ{2) ■ zi "J^'p" 

- + (£!y), + (Zr q]q. 

\ ka| Q V pal/ a / 



U 13 



(20) 



I •2a I I -^■a I 



(21) 

For (20) we write 



\ztf = Q'^lBRf + C^\Za\'^ + 2Q2cBi? • Za 



Now 



which means 



We write 



(20) = NICE - ^' I I Qa, 



(20) + (21) = NICE - - 2QQc«Si? • 



Zat — {Zoit 'Zee) I 12 {^at ' •^■q. ) i |2 
only depends on t 

= 5(t) Za + ((Q'Si?)„-4 + CZaa-4)^ 
See (IV.21) 

= B{t)Za+ D{a) ^-^Zaa-Z^r^. 
as in (IV.45) 
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Writing zt = Q^BR + cza we compute 

Zat ■ Zt = Q^BR ■ ZaB{t) + DQ^BR ■ 



NICE 



NICE because D is nice 



- T^Zaa ■ Z^Q^BR ■ ^ + cB{t)\Za\^ . 
\Za\ I Za I V 

NICE 

To simplify we write 

Zc,t ■ Zt = NICE - j^Zaa " z^Q^BR ■ j^. 

I Za I I Za I 

Setting the above formula in the expression of (20)+(21) allows us to find 

(20) + (21) = NICE . 

This yields 

(fiat =NICE - I -(fiaa 

\Za\ 

- Q2 ( (BRt + J^BRa) + (zat + J^Zaa 

\\ \Za\ ) ^ZaV V \Za\ 



+ V5P2"'(^))-4 

- Q' ( \BRV + ^ + 2c^^ - ^ ) VQiz) . 4 """ ■ 



+ 2c^^^ - ^] VO(z) ■ z^'-^^ 



+ {Quj + 2QBR ■ Za)VQiz) ■ z^ ' ^' 



\z P ■ 

We now complete the formula for a in (IV. 2) to find 

<^at =NICE - 1 -Apaa -Q CT- 



, Za I I ■^a I 

3 



+ Q 



2 



''a 



(22) 

+ ( -l^i^l^ - ^ - + ^) VQ(.) . 4^ 



(23) 

+ {Quj + 2QBR ■ Za)VQ{z) ■ z^^f^ 



(24) 

Expanding 
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we find 



(22) + (23) = Q\^.+BR.z.^^- 2c^ - VQ(.) • 
Writing 

Zat ■ = {Cf'BR)a ■ + CZaa ' Z^ 

we obtain that 



(24) = {Qco + 2QBR ■ Za) VQ{z) ■ z^ 



\z |3 



+ {Qu + 2QSi? • Za) VQ{z) ■ z^c-^f^ 



Thus 

(22) + (23) + (24) = g3 + BR • Vg(.) • 

.(Q25i?),.2;^ 



+ (Qw + 2QBi? • Za) VQ(z) • z^ 



I ^OL I 



QVQ(.) • 4 (c + 1BR . Za) ( + ^^''f'^j" • "° 



2lz P Iz |3 



QVQ{z) ■z^{oj + 2BR ■ z^) T^D{a) 



= NICE. 

Finally, we obtain 

^at = NICE - ^^aa " Q^^^f^. 

\Za\ I ■2a I 

□ 

Corollary IV. 13 // we disregard the condition on the H^~^ norm for the definition of the 
NICE terms, imposing only the first condition, then 

= NICE-Q^a- 



IV.D.ll Higher order derivatives of a 

In this section we deal with the highest order derivative of the R-T function. We show that 
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Lemma IV. 14 Let z{a,t) and (jj{a,t) he a solution of (1. 19-1. 20). Then, the following 
identity holds: 

di-\Q^a) = ANN+ \za\H{dt'^t) + ^H{dl^) (IV.48) 

where ANN satisfies 

\\ANN\\l2 < CEP{t) (IV.49) 
for k > 4, where C and p are constants that depend only on k. 

Proof: We show the proof for k = A. Prom now on, if a term / satisfies 

Wfh. < CE^{t) 

we say that this term becomes part of ANN. By abuse of notation we will denote / by 
ANN. We recall 



gV =Q2 ( BRt + ^BR^ ■ + (zat + ^ 



+ 



BR^ 



2 1 Zrv I ^ 



VQ-z^ +Q'VgP^\z)-zt. 



this term is in so its third derivative is in ANN 



We write 



Zoit \ T Zci 



Zq\ 



2 1 Zct I ^ 

2 1 P 

2 1 Z(x I ^ 
2 1 Zqi I ^ 

Above we use (IV.45) and (IV.46) to find 



this term is also in 



{Q BR^oi ■ Zq, + CZaa ' ^~ i i ^oia ' Z^ 

Zn 



iQ^BR)a ■ z^ + 



D(a). 





Za 1 






2 1 Zct 1 ^ 



Zat + 



.Z^ = AN+^//(V.a), 



'^\Za\^ \ \Za_ 

where AN is as in (IV. 47). The remaining terms in Q^a are 

L = Q^BRt ■ 4 + j^BRa ■ z^. 

\Za\ 

We take 3 derivatives and consider the most dangerous characters: 

dl{L) = Ml + M2 + Ms + ANN, 



(IV.50) 
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where 



Ml = Q^BR{z, dltJt) ■ 4 + ?^BR{z, 5^a;) • z^, 



J —IT \^ Z I 

\zJ 27r / ^ \z - 2;'I2 



M3 



■a\T^ J —TT \Z Z I 



Here we point out that in order to deal with BRt in the less singular terms we proceed using 
estimate (IV. 30). In M2 we find 

M, = ^Aidlzt . z„) + ^A(9^ • ze.) + ANN. 
For the second term we use the usual trick 

d^Z ■ Za = ^9^Z ■ Zaa- 

For the first term we recall that 

\Zaf = A{t) Za ■ Zat = ^^'(*) =^ (^a ' ^at)a = 

^ Zaa ' Zat ~l~ Za ' Zaat — ^ Zaaa ' Zat ~\~ "^Zaa ' Zaat ~^ Za ' Zaaat ~ 
=^ Za ■ Zaaat — "^Zaa ' Zaat Zaaa ' Zat- 

This allows us to control M2. For M3 we find 

M3 = -f^A{za ■ dlzt) - f^Aiza ■ diz) + ANN 
I Za I I Za I 

SO it can be estimated as M2. There remains Mi. Using that {z — z')-^ ■ z^ = {z — z') ■ Za 
we find 

Ml = ^Hidlut) + ^^Hidtco) + ANN. (IV.51) 

We compute 



2 /i. 

+ H{dl{\za\c)t) + ANN 

+ H{dldt{-{Q^BR)a ■ Za)) + ANN. (IV.52) 
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We compute the most singular term in 

q2 {dlzt-dlz[ 



•^TT \Z Z \ 

71 I'Z Z \ 



extra cancelation in (2— z')-'--2(j=(z— z'— z„/3)-'--2(j 



27r 



extra cancelation as above 



This shows that 

dldt{-{Q''BR)o^ ■ za) = -^^k{dlzi ■ z^) + ANN. 

That gives 



dldti-iQ^BR)^ ■ za) = -A ( S^diz^ ■za]+ ANN, 



2\Za 



which imphes 



H{dldt{-{Q''BR)^ ■ z^)) = (^^a^t^ • ^a) + ANN = (dlzt ■ z^) + ANN. 

Plugging the above formula in (IV.52) we find that 

^H{d',ut) = IzalHidl^t) - (dlzt ■ z£) + ANN 

= K\H{dl^,) - {dliQ^BR) . z£) - (c5^ • z£) 

+ ANN. 

As we did before, we expand da{d^{Q^BR) ■ z^) to find 

d^idUQ'BR) . z^) =2QVQ{z) ■ dtzBR ■ ^ + ^PV J ^^^^0^^u;'d/3 

Therefore, we can use (IV.36),(IV.37) and (IV.37) to show that the most dangerous term is 
given by Q'^\H{d'^u). It implies 

d^{dl{Q''BR) . 4) = Q^\H{diu) + ANN 
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and therefore 



Midgut) = \z^\H{dl^,) - ^^2^H{d^a^) - (^^ • 'a) + ANN. 



We use the above formula and expand (p to find 

Ml = ^Hidlu;^) + + ANN 

= \za\H{di^t) - ^cHid^oo) - ^^cd. {d'^z ■ z^) + ANN 

= \z^\H{dl^,) - cW\H [dt (gj)) - ^ca. [diz . 4) + ANN 

= \Za\H{dlipt) - c\Za\H {d^^ip) - c\Za\H{d^^{c\Za\)) - ^^^^a (d^ ■ Z^) + ANN. 

We will show that 

-c\z^\H{dt{c\za\)) - ^(^da {dtz ■ 4) = ANN. (IV.53) 

It yields 

^Hidlcot) + ^^H{diu) = \zamdl<pt) - c\za\H [dtv] + ANN 
that together with (IV.50) allows us to obtain (IV. 48). We have 



-c\z, 



.\Hidi{c\z^\)) - ^Cda (d^Z . Z^) = -cH{d'MZa\')) - ^cS^ • 4) 

= cHidliiQ'BR)^ . z^)) - (dtz ■ z^) 

We repeat the calculation for dealing with the most dangerous terms in 
dlUQ'^BR)^ ■Za) = A (dtz^ ■ za^) + ANN. 



We recognized as before terms in ANN using that {z — z')-^ ■ Za gives an extra cancellation. 
We find that 

cH{dU{Q'BR)^ . z^)) - ^c9« {dU ■ z^) 

= ch(a( dU^ ■ Zo.^\\ - S^cdo. i dU ■ z}-'] + ANN 



"'^OU 12 / / oU |2 ^ V ^ 

^\^OL I / / ^\^OL I ^ 



= -Cda d^Z-^ ■ Zn 



Using that d^z-^ ■ Za = —d^z ■ z^ we are done proving (IV.53). □ 
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IV.D.12 Energy estimates for 

In this section we prove the following result. 

Lemma IV. 15 Let z(a,t) and u{a,t) be a solution of (1. 19-1. 20). Then, the following a 
priori estimate holds: 



d 



m 



(t) < -Sit) + CEP{t) 



for k > A, where C and p are constants that depend only on k. 



(IV.54) 



Proof: We shall present the details in the case = 4, leaving the other cases to the reader. 
Using the estimates obtained before one has 



\ip\\Ut) < CEP{t). 



Developing the derivative using Lemma IV. 12, we get that: 



d^ 
di 



=h + h + h, 



(IV.55) 



where 



h = 2 A{d'lip)dt,{mCE)da, h = -2 A{d-i^)di{^^^aa)da, 

11 Jj \Za\ 



cr— 



)da. 



We use (IV.41) to control Ii. The most singular term in I2 is the one given by 
-2 / j^A{dl<p)dt,^<pda =2 I J-AV2(93^) Lk^l^dl^) - K^'^^di^) 

JT \Za\ JT pal '- 

+ / ^^dMK^/\dl^)\''da. 

Jt \Za\ 

Using the commutator estimate 

\\gA'/\d^f) - A'/H9daf)\\L^ < \\g\\c4f\\H^r^ 
we can bound l2- In we split further considering the most singular terms 



da 



(IV.56) 



Ji 



Jt 



J2 = -2 / A{dl^)Q'a^ 



2 L/q,^ 



da, 
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J3 = -2 / A{dl^)dl{Q'af-^^da. 



The term Ji can be estimated as before. Recalling (IV. 35) we see that J2 = —S{t). It remains 
to control J3 in order to find (IV. 54). 
We decompose J3 = Ki + K2 where 



and 



K2 = 2[ H{dl^)dUQ'af-^f^da. 
Inequality (IV.31) for A; = 4 allows us to obtain 

Ki < CEP{t). 

To finish the proof we use formula (IV.48) for A; = 4 to find K2 = Li + L2 + Ls where 

JT \Za\ 

L2 = 2 [ H{dl^)\z^\H{dl<p,f-^^da, 

JT l^a] 

L3 = 2 / H{dl^)<pH{dtvf-^^da. 

JT pal 

The term Li can be easily estimated using (IV.49). For L2 we substitute the expression 
(IV.54) for d^cpt to get L2 = Mi + M2 + M3: 

Mi = 2 / H{dl^)\za\H{dl{mCE))^^^^^da. 

JT \^a\ 

M2 = -2[ Hidl^)\za\Hidli^)f-f^da. 
M3 = -2 / H{dl^)\z^\Hidl{Q'a'-f^)f-^^da. 

JT \^a\ \Za\ 

By equation (IV.41), Mi is bounded. M2 is bounded knowing that we have room for half 
derivative in the term which is not the third factor. Finally we can bound M3 in virtue of 
Lemma IV. 8. To finish, in L3 we integrate by parts to find 

L3 = - / \H{dl^)\^daif^^^^^)da<CEP{t) 

JT \^a\ 

using Sobolev embedding. □ 
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IV.D.13 Energy estimates for . " + V , 

Lemma IV. 16 Let z(a,t) and uj(a,t) be a solution of (1. 19-1.20). Then, the following a 
priori estimate holds: 

^ ( ^ + y < CEP{t) (IV.57) 

dt\m{Q^a){t) ^m{q^){t)/ ~ ^ ^ 

for k > 4, where C and p are constants that depend only on k. 

Proof: Inequalities (IV.27) and (IV.32) show that (Q^a) G C^([0, T] x [-tt, tt]) for some T and 
therefore m{Q'^a){t) is a Lipschitz function differentiable almost everywhere by Rademacher's 
theorem. Let 

m(QV)(t)= min {Q^a){a,t) = {Q^a){at,t). 

7r,7r] 

We can calculate the derivative of m{Q'^a){t), to obtain 

(m(QV)y(0 = (QV)t(at,t) 
for almost every t. Then it follows that: 

d f 1 \ iQ'^cr)tiat,t) 



dt V"i(QV) J ' ' (m(QV))2(t) 
almost everywhere. By using the previous a priori estimates for the L°° bounds, we get to 

d ( W? \ ^ CEP{t). 



dt \m{Q'^a) 

On the other hand, wc can apply the same argument to — . . Denoting again by at the 

m{q^)[t) 

point where the minimum is attained we have that: 

df 1 \ 1^^^ _ zt{at,t) ■ {z{at,t) - q^) 



dt \m{q^) J ^ ' {m{q^))^{t) 
which again can be easily bounded and we get (IV.57), as desired. 

□ 

IV.E Proof of short-time existence (Theorem 1.7) 

To conclude the proof of the local existence, we shall use the previous a priori estimates. We 
now introduce a regularized version of the evolution equation which is well-posed for short 
time independently of the sign condition on a{a, t) at t = 0. But for a{a, 0) > 0, we shall 
find a time of existence uniformly in the regularization, allowing us to take the limit. 
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Now, let z^'^'^{a,t) be a solution of the following system (compare with (IV.20)): 

(IV.58) 



' Q'^iz^'^'i^y ^ \\^az^'^'^'{a,t)\ I 

- 2Q{z^'^''')Q{z^'^'^')a\BR{z^'^''',0J^'^'^')\'^ - da 



2BR{z^'^''',uj''^''')f zl;^^" 



+ 



r {Qiz^'^'^'fBRii 

J —TT 



7rQ(z^'^'A')2 



d/3 



2\^az'^^^^'{a,t)\ 



2a 



+ 



2e 



(IV.59) 

'^'^(a, 0) = zo{a) and a;^'''''^(a, 0) = oJo{a) for £ > 0, (5 > 0, /x > 0, 05 and 0^ even mollifiers, 



and 



27r 



/a 
-71 



05 * 05 * (5^(g2(^^'^''^)(^)i?i?(z^''^'^a;^'^''^))(/3))d/3, 



2\daZ^'^^l'\ 



B^'^'i'it) = ^ r if"5l^r/"'Ii9 • da{Q'^{z''^^'')BR{z''^^i',uj''^''')){a,t)da, 



* 05 * 



27r y_^|5a2^'^'/^(a,t)|2 
+ 7r r 9^z^'^'^(/3) 



2vr |a3Z--.'5'^(/3)| 
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We start proving the following lemma: 

Lemma IV.17 Let z'''^'''{a,t) G H'^iT), u^'^'i'{a,t) G H'^iJ), ip^'^^''{a,t) G H^{T). Then 

Lo'^^'^'{a,t) G H^{T). 

Proof: We can write ul^'^'^^ as: 

, ,£,<5,/i _ "^IdaZ^'^'f^l f e,S,iJ. , r-eAA 

Taking three derivatives yields 
dlu^'^'^ = SAFE + ^l^a^C^'^-- 



SAFE - * 0. * (^1^^ • (Q^(.-^"^)a3si?(.-^-,a.-^-))) 

SAFE - f^^. * 0. * • 53.^.^-))] 



where SAFE means bounded in L^. Using the representation 



BR{z'^^''',dlu;'^^^'')= SAFE + ^ ^,^"^1 g^l Hidlu''^'") 



we get that 



a^o; = SAFE - Q^^<Ps *<Ps*[Q{z' '%^Q^;^Hidg^ ' \d^\) 

=0 

and we are done. We should remark that the lemma holds independently of (5, /x and e. □ 

We define a distance between data {z,ijj) and {z,uj) by taking 

d{{z, Lo), {z, Ld)) = \\z — z\\fj4, + — ^^iWh^ + \\<f — ^\\h^ 

where (p and ip arise from {z,oj) and {z,lo) respectively by (IV. 18). Let XX denote the 
resulting metric space. The proof of Lemma IV.17 gives also the following 

Corollary IV.18 The map {z,oj) ^ uj is Lips chitz from any ball in XX into H^{T). 

We note that throughout this section we will repeatedly use the following commutator 
estimate for convolutions: 

Us*{daf9)-g<Ps*{daf)\\L^<C\\dag\\L'^\\f\\L^, (IV.60) 
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where the constant C is independent of (5, / and g. We can now operate to get the following 
expression for ip^'^'^^: 

~)2( e,5,ij.\ p 



2\^az^^^^^'{a,t)\ 



The RHS of the evolution equations for z^'^'^ and ip^^^'^ are Lipschitz in the spaces i?^(T) 
and H^'^^i^) since they are mollified. For the case of u^'^^^ (Lipschitz in the space i7^(T)) 
we use that for 5 small enough (j)^ * </*<5 is close to the identity and the a priori bounds. In all 
of the cases we have taken advantage of Lemma IV. 17. Therefore we can solve (IV.58-IV.59) 
for short time, thanks to Picard's theorem. 

Now, we can perform energy estimates as in the a priori case to get uniform bounds in /x 
and we can let jj, go to zero. The energy estimates that we can get are the following: 



dt 



it). 



We should note that for the new system without the (f)^ mollifier, the length of the tangent 
vector |9q,z^''^| is now constant in space and depends only on time. Lemma IV.17 still applies 
and we can still perform energy estimates as in the a priori case. The only difference relics 
on the fact that we should have to move the mollificrs and apply the estimate (IV. 60). We 
should also remark that because of the dissipative term eAtp^'^ it is enough to use the following 
estimate 
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and hence require only that dt^p^'^ G H^{T) (instead of the H''^+2{T) that was required 
before) except for the transport term that can be estimated as in subsection IV.D.12. The es- 
timations are performed following exactly the same steps of subsection IV.D. More precisely, 
we can get the following energy estimates: 



di 



f 1 A 



< c{e) \\z^''\\%. + ||^(.^'^)||ioo + \\co^'%. + y^'^W'. + E z::eJ7:r. (*)• 



^ if ' ] 

^^^^ m^'^il^) J 



3 



Under these conditions, we can let 5 go to zero. 

Finally, let z^{a,t) be a solution of the following system (compare with (IV.20)): 

zl{a,t) = Q^{z'){a,t)BR{z' ,u^){a,t) + c^{a,t)daZ%Oi,t), (IV.61) 



Qiz-r. 

+ nQ(z^^ Jj^^' ^ ' ^^'^ ■ W) ^ ^ ^ 

z^{a, 0) = zo{a) and oj^{a, 0) = ujo{(x) for e > 0, where 
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Proceeding as in section IV. C (compare with equation (IV.43)) we find 

daipt = -B {t)ip^ - + da T^— JT7^7-7y(</=' ) 



2\do,Z^\ ''\\daZ^\Q{z 

- Q{z'fdtBR{z',u') • d^z' ^^^f • ^^l^' - dt{\do^z'\B^{t)) 

+ daiQ{zn'BR{z',iv')) ■ g^^^ ^^y \^^'\ 2daiQ{z')dtQ{z')BRiz',u')) ■ ^ 



da ( 2c^BR{z^,u^) . ^^Qiz^)daQiz^)^ - d^ 



We also define (compare with equation (IV.2)) 



(IV.63) 



a' = {dtBR{z',u') + ^^daBR{z',uj')) • d^z' + l^^(5«zf + ^fl^') ' ^a^' 



+ Q(^^) 



2 acwZ^ r 



' (VQ(z^)) • +5(VP2-^(z^)) • 



Remark IV. 19 The system (IV.61-IV.62) is analogous to the system considered in [13, 
Section 8]. We point out an unfortunate typographical error in that section; the Laplacian 
should have been written as the square root of the Laplacian. 

For this £-system (IV.61-IV.62) we now know that there is local-existence for initial 
data satisfying T{zo){a, P) < oo even if a^{a,0) does not have the proper sign. In the 
following we shall show briefly how to obtain a solution of the regularized system with 
z' e C{[0,T'],H''),(p^ G C([0,T^],i7^-5),a;= G C{[0,T^], H''-^) for k > 4. 

The next step is to integrate the system during a time T independent of e. We will show 
that for this system we have 

j^E{t) < CEP{t), (IV.64) 

where E{t) is given by the analogous formula (IV.23) for the £-system, and C and p are 
constants independent of e. 

In the following we shall see what is the impact of the e system on the a priori estimates 
and check that there is no practical impact for sufficiently small e. To do that, we will show 
the corresponding uniform estimates for k = 4 and leave to the reader the remaining easier 
cases. Let us consider the one corresponding to I3 in section IV.D.12, we have 

II = -2 j^,A{dl^%a))dliQ\zna^dlz^ ■ d^z^)da. 
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Proceeding in the same way as before, we can perform the same spHttings and get uniform 
bounds such that if = — + M| + "bounded terms" where corresponds to S in (IV. 35), 

I "bounded terms" | < CEP{t), 

and 

M| = -2e r ^^^f H{dl^^)H{Kdl^-)da. 

J— TV \'-'a^ I 

Then we can write M| as follows 

Ml = -2e £^Ah(^^^-^H{dl^^))A^Hd'^^')da, 
and therefore, for small e 

Ml < ||Ai(9^(^^|||2 + "bounded terms", 

which gives 

f^E{t)<CEP{t)-'-\\A{dl^)\\l,. 

This finally shows (IV.64) and therefore 

E{t) < {Ct{l-p)+E^-P{0)Y/^^-P\ 

Now we are in position to extend the time of existence so long as the above estimate works 
and obtain a time T dependent only on the initial data (arc-chord, Rayleigh-Taylor, distance 
to the points q^, . . . , q^, and Sobolcv norms of z,u}, and 99). We can let e tend to 0, and get 
a solution of the original system. This concludes the proof. 
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